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In  a  previous  project  a  class  of  RF  intrusion  sensor  systems  was  analyzed  and  numerical 
results  were  obtained  for  various  intrusion  situations.  The  system  consisted  of  a  ported 
coaxial  cable  laid  circularly  around  an  area  to  be  protected  and  an  antenna  near  the 
center  of  the  circle.  The  presence  of  an  intrudeV  (electromagnetic  scatterer)  near  the 
cable  perturbs  the  field  and  thereby  changes  the  signal  at  the  antenna.  The  analysis  was 
based  on  the  plane-wave  spectral  representation  of  fields  and  the  assumed  scatterer  was  a 
uniform  spheroid  with  constitutive  parameters  designed  to  simulate  those  of  a  human  frame. 
Emphasis  was  on  polarization  effects. 

Although  the  previous  analysis  constituted  an  attempt  to  model  this  system  quite  accurately, 
the  price  for  this  accuracy  was  extremely  large  computer  time.  In  order  to  reduce  computer 
time,  approximations  were  made  along  the  way  that  reduced  the  overall  accuracy.  The  present 
report  discusses  a  follow-on  to  the  above  mentioned  project  whose  objective  was  to  produce 
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a  simpler,  less  computer-time-intensive  model  and  associated  software  that  ran  be 
readily  used  to  study  the  effects  of  parameter  changes  on  system  performance. 

Steps  were  taken  toward  that  objective,  as  follows: 

1'  The  software  was  thoroughly  checked  and  redesigned  for  greater  efficiency 
when  that  step  was  warranted. 

'21  The  model  for  the  cable-generated  fields  was  simplified. 

,3)  While  many  of  the  features  of  the  plane-wave  spectral  representation  of 
fields  •■•ere  retained,  that  approach  was  compromised  by  approximating  the  cable- 
generated  fields,  within  the  scatterer  voluae,  as  those  of  a  single  plane  wave 
propagating  in  the  direction  of  power  flow  of  the  true  cable  generated  fields.., 

Th“  result  of  these  modifications  is  a  new  computer  program  with  much-reduced 
running  time  which  produces  numerical  results  chat  are  much  more  consistent  and 
interpretable  than  chose  produced  in  the  previous  project. 
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\ .  Introduction 

The  work  described  in  this  report  is  a  continuation  of  research  done 
nn  RAT  Postdoctoral  Contract  F30602-7ft-C-0Kt2,  completed  on  ‘November  1,  IP81. 

The  Tinal  Report  on  that  contract,  entitled  "Analysis  of  the  Polarization 
Vpendence  of  the  Interaction  Between  Human  Frame  Targets  and  Radio  Frequency 

„  .-13 

Sensor  Melds,  is  desit?nated  as  Reference  1  in  the  present  report. 

■^e  previous  project,  referred  to  above,  involved  nalysis  of  a  particular 

tvpe  of  PF  intrusion  sensor  system.  The  system  consist;  of  a  ported  coaxial 

cable  laid  circularly  on  the  ground  surface  and  a  verti'^al  receiving  antenna 

at  or  near  r'^e  center  of  the  circle  tuned  to  the  freque'cy  of  the  wave  propa- 

patina  in  the  cable.  (See  Figure  1.1).  ^n  intruder  stepping  or  crawling 

icross  ’.he  c.'hl®  in  order  to  enter  the  circle  scatters  '•ome  of  the  electro- 

nagretic  energy  Impinging  on  hix  or  her  from  the  cable  into  the  direction 
2  4-6  S 

of  the  incenn.i,  changing  the  signal  received  by  the  antenna  and  thereby  detecting 
the  presence  of  Che  intruder. 

although  the  concept  of  this  device  is  very  simple,  the  electromagnetic 
Choorv  problems  encountered  in  attempting  to  model  its  operation  analytically 
are  extremely  complicated.  In  the  original  project,  as  renorted  in  Reference 
1,  Che  scheme  upon  which  the  an.olysis  was  based  was  the  plane-wave  spectrum 
reprps^ntat  ion  of  an  electromagnet  ic  wave,  '^eelectrir  n(i  magnet  ic  fields  fror  the  rahlr 

V  ere  calculated  based  nn  analysis  of  the  propagation  al  ng  the  cable  and  the  use  of  the 

«  6 

"FtntLon-rhu"  (or  "Kirchoff  Huyghens’*)  integral  formui.;  -hich  gives  the 
fi-^ld  comnonents  in  space  based  on  a  knowledge  of  the  fields  at  the  cable 
Si^ts.  ’’his  was  done  as  if  the  cable  '-'cre  in  free  space.  Then  the  plane- 
’.avA  spoctrnl  representations  of  each  the  field  comp'-nents  was  calculated. 

"p’.nnc  '  ivo"  was  then  taken  through  the  ground-reflection  process. 


Figure  1.1 

Intrusion  Sensor  ^vstem  Configuration 


( direct -:)lur.-«rf>'jnf!-reflected)  spectral  field  in  the  absence  of 

-■•.n  intruder  wrs  calculated.  The  intruder  was  "lodelled  as  an  electromagnetic 

scaiterer.  Tlie  t'^  tal  spectral  field  at  *^he  intruder’s  location  vas  ronsidered 

as  the  field  o:  a  plane  wave  ineiJent  on  the  scattercr,  which  was  modelle! 

IS  a  uniform  spheroid  with  constitutive  rnraneters  chosen  from  a  composite 

'f  those  of  the  v  rious  constituents  of  the  human  body,  \  pronrar.  origirnllv 
1,1,  45-57^50, 5?, 57, 

developed  }  y  Pete’  Rarber  at  the  I'niveisity  of  I'tah  was  used  to  perform  the 


scnttcrine  analys  s.  This  progran  was  inserted  into  our  so‘‘tware  as  i  siibrojtir 


Tne  scatter^i  wnvp  fields  propagatirn^  directly  toward  rhe  artenna  and 
those  re‘‘lected  from  the  ground  and  then  propagating  toward  the  antenna  were 


'"O’^puted . 

The  program  developed  for  this  problem  computes  the  components  of  tie 
plnne  wave  spectrum  of  four  contrihiitions  to  the  electric  field.  They  are: 

Ca"!  the  field  from  the  cable  in  the  absence  of  the  ground  (i.e.,  in 
infinite  free  space) 

''b)  the  ground-ref lected  cable-generated  field 

(c)  The  field  scattered  by  the  intruding  object  directly  into  the  antenna 

^d)  the  ground-reflected  scattered  field  in  the  direction  of  the 
antenna 

•'ach  of  the  constituents  (a),  (b).  (cl  and  (d)  is  a  function  of  the  normalized 
wave  propagation  vector  The  final  step  in  the  computation  v\ns  the  invervSe 
f'^urier  transformation  of  the  spectral  electric  field  components  in  order 
CO  produce  the  field  components  at  the  antenna  site  as  a  function  of  the 
nosirion  of  the  antenna.  The  end  result  of  the  computations  wrs:  (1)  a 
of  field  components  consisting  of  the  superposition  of  contributions 
and  fh\  i.e.,  the  field  conpon-.nts  in  the  absence  of  the  scattercr; 

'21  tho  sunerposit  ion  of  contributions  (c)  and  (d),  i,e.»  the  field  conpon- 
'•'rir-s.  luo  to  scatterer;  and  ^11  the  superposition  o*"  concrihut  ions  (nl, 

'  ),  'c'  and  ^d),  i.e,,  the  cornoneats  of  the  total  electric  field  seen  at 
tho  antenna  fice. 


Ihr  purpose  of  the  fnl low-on  project  was  to  produce  an  improved  analytical 
“odol  in  improved  computer  program  that  would  be  more  useful  in  studving 
svstem  r^rformnnee.  The  problems  with  the  original  analvticnl  model  and 


assoc  I'lt'’ '  ‘^oftKirc  ’>er'': 


1 . 


I'  U  Ip.'hs*  -5"\r"3Tnl  and  \n  ?one  c.-iros  qucFt  ionnblr .  This  co'*',proni  sos 
ncrurncv  of  the  conputat  ions , 

f  ‘  "^0  the  dif  f c*jl  t  ies  ^encountered  on  the  original  pro'cct, 

*  o-  w«_‘ro  'ir  “’v  iue  to  the  attempt  to  put  the  analysis  of  this  s‘-3tc~' 

a  sr'l'.  ‘  fteorrt  ical  foundation.  e  were  very  careful  about  the 
ssannt  io:v  invo';el  anti  trieJ  to  ■aa'-:e  the  nodel  as  realistic  as  possible, 
•'.ever,  ’.ero  forced  into  trade af between  the  realism  of  the  Todel 

'  i  t^o  Ijmitati-es  o‘  comnater  tioe.  Although  the  equations  proqrar.re! 

•  r  •“•‘c'.'.ae  rannutation  may  hav''  contained  a  great  deal  of  realism, 

^  to  save  running  tire  (e.g.,  insufficient  points  iJi  ca  murerical 

inr--' prat ’.( n)  might  in  <ore  case*^  produce  a  result  lo?s  accurate  than 
'-■-’ill  navo  boon  o'^toined  v.it‘'  a  le<^s  rigorous  hut  accurately  corrutablr 

“’“'isurea  tako**  on  the  neu  '^roiect  to  improve  the  analysis  tnd  corpu- 
aMons  are  'iscussoc  in  the  present  report,  fn  Section  2,  t'-'o  calculation 
■  tbe  from  the  c<able  slot^^  as  if  the  cable  were  in  free  speoce  are 

Irscrihe;.  In  tao  oariier  project,  these  caicuiatioms  were  made  for  erbitrarv 
i''-les  prop<a:afing  along  t!ie  cable,  i.e,,  TF!,  TM  or  TF*(  .modes  (?ef.  I,  Appendix 
'  bilo  that  model  was  very  general  and  hence  wouhl  cover  many  possible 
■pses,  It  -..'ns  also  '•iimbcrsonie  end  required  j  large  anounr.  of  so^twnre. 
as  doc  i  led  t''  ‘'ormulnte  the  field  s  on  the  coble  slots  using  a  simpler  model 
ose'^  f'U  orooncotion  of  a  ”T  I  mode  jlong  the  cable.  This  yields  a  simpler, 

•’ere  interpr-^table  sc*  of  '*ornuIas,  requires  considerable  less  softi.iro 

md  vjf!  fc  :'^ntlv  rrnijstir  r'or  the  purposes  of  thi*^  analysis.  Also,  it 

.Mcilv  he  gcnerr.l  i  7ei’  if  nece-^'^iorv  and  if  data  ver^  .available  to  allow 
i.S'-’'  rmort  of  pararet^'r  values  *or  !  -an't/or  T'  modes. 


-urre  * 


In  section  3,  the  '■emainino  steo*?  ir.  calculation  of  tho  rcceiv^fi 
at  the  vitprna  are  described.  In  effect*  Section  3  contains  brief  su“'^Dries 
or'  th"  sections  in  Reference  1.  Section  3.1  briei'lv  surrj^ar i.Tes  the  covrro;? 
of  til'?  plane  V'ave  soectral  representation  of  fields  that  was  presented  ro 
?reat  detail  in  Sections  3  and  4  and  Appendices  il  and  III  of  Fef.  1.  The 
enpnasis  in  'Section  3.1  is  on  the  najor  results  that  vere  derived  in  llefercncc 
\.  So“re  nev  investigations  on  the  use  of  this  approach  are  introduced  in 
Section  3.1  and  further  ‘etailed  in  Appendix  II. 

Section  3,2  contains  a  brief  sunmar>  of  the  key  points  on  the  grourri 
reflected  wave  that  were  covered  in  ereat  detail  in  Section  S  of  1. 

‘o  significant  changes  in  the  analysis  or  the  propraras  were  nade  on  the  pre¬ 
sent  project, 

'Sections  3.3,  3.4  and  3.>  contain  brief  sumraries  of  the  'oateriol  prerenter 
in  detail  in  Sections  S,  7  end  8  respectively  of  “^ef.  1,  This  naterial  concernr 
the  coorrinato  transfor"\ations  to  allov  the  superposition  of  direct  and  '■  ronnd- 
refloctel  rcble-pcnerated  fields  to  be  rodelled  as  the  field  incident  on 
the  sentterer  (Section  3.3),  the  action  of  the  scattering  proprar^  itrol:' 

(■■ection  3,4)  and  the  transf ornat icp  of  the  scati'^red  field?  fron  scotcercr- 
centered  coordinates  to  earth  coordinates  <'!^eciion  3,3).  Soclien  3.h  is 
surnarv  .if  prour.d-ref lections  of  V:c  .scattered  fields,  prer-ente  !  in  det.M  I 
1  "cr  t  i  }  of  ‘e : .  I . 

Section  3.7  is  a  sunnary  of  the  process  of  sunminit  all  field  contri  hut  ions 
and  taking  the  inverse  Fourier  trans''or':  in  order  to  convert  ^ro"  tlie  plano- 
’•ave  spectrur  of  the  field  conponents  to  the  actual  field  connonents  as 
functions  of  «pace  coordinates. 

'ection  4  rontcins  a  detailed  .nalvsis  of  the  new  aporoach  t nt  vns 
finally  n::ontr‘i'  the  1  jfer  ''roj-’ct.  ‘his  c.pproarh  hn.r  r'^''l-nced  t*'':  ol.inc-.'cvo 


.V.-* 


•  ’  •  *  W.  •  •  "v ' 


’.v.-Sv'.v 


!=*iortr«n  .ipproach  but  ba''-  r<>tain»>d  nost  of  its  features.*  It  involves  the 
apnroxination  of  the  field's  ioci^ont  on  the  scatterer  as  those  of  a  single 
nine  wave  'not  a  spectruni  of  pl^ne  waves)  whose  direction  of  propagation 
' that  of  the  Povntin'j  vector  (direction  of  newer  flow'  jt  the  center  of 
t'^e  sr  at  r  er'^r .  This  approach  c  i rr  .nrvents  the  proble’r^  of  performin':  the  two- 
*  i-^eosionil  ’averse  Pourier  transformation  and  hence  saves  an  enormous  anount 
of  corputor  time,  allowint:  many  more  parameter  regimes  to  he  studied  ’./ithir. 
a  'lives  timo.  Ir  view  of  the  degree  of  approximation  required  to  implement 
tho  co'*‘n'U''’t  ions  for  the  plane-wavc-wave  spectrum  approach  f^ee  Appendix 
n  fer  detailed  discussion  of  these  problems)  the  new  approach  offers  many 
Tivantages  in  addition  to  that  of  increased  computational  speed.  If  the 
rel-it ionsh ips  arising  in  the  plane-wave  spectrum  approach  could  be  comnuted 
TCC'jrately  without  a  prohibitive  expenditure  of  computer  tine,  that  would 
be  1  sjoorior  methodology.  However,  the  required  compromises  in  computational 
ac'‘’iracv  night  negate  the  advantages  of  the  theoretical  model  Itself.  The 
nmroarh  involves  an  approximat ion  at  the  initial  formulation  stage  which 
is  phvsic.illv  satisfying,  l.e.,  the  idea  that  the  **infident  plane  \iave"  that 
anproxi-intos  the  actual  wave  in  the  scatterer  region  is  a  linesrly  polarized 
rra'^sverso  electromagnetic  v/ave  traveling  in  the  direction  of  the  power  flow 
fh‘>  fields  emanating  ^rom  the  cable  at  the  scatterer's  center, 

■^I’.e  new  approach  orodiiccs  results  that  appear  to  he  very  reasonable 
ir.-:  ph'.sirallv  internretnble.  Since  comuter  tine  is  not  so  critical  with 
this  approach,  it  is  possible  to  study  many  variations  of  parameters  within 
n  reasonable  expenditure  of  tine.  All  of  this  is  discussed  in  detail  in 
Section  5  of  this  report. 

Ihe  normalised  wave  nropagation  vector  which  was  used  extensively  in  the 

apnroarh,  is  still  retained  in  the  new  approach.  It  is  used  in  nodcllin^ 
of  -.rnftnring  and  ground  reflection.  Ml  stages  of  the  analysis  reported 
.n  -ection  3  remain  ossontially  intact,  with  only  minor  modifications  required. 


..  4 

• '.N 


2.  Calculation  of  rielHs  fron  Cable  Slots  ir  "roe  Soacc? 


Tho  an'^rnach  used  in  the  previou**  project  (Ref.  1,  Sections  3  and  4, 

and  Apncmiicrs  ITT  and  IV)  has  been  considerable  sinplified.  In  the  revise'^ 

nnnroacb,  the  ronponencs  of  Che  field  frorn  the  slots  is  calctilated  directly 

in  the  spatial  do'iai''  rather  than  in  the  spectral  do-nain  (Ref.  1,  Appendices 
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II  nrd  in)  as  in  the  previous  «orl«'.  Av^ain,  the  Stratton-Thu  (or  Airchoff- 
•hivphrns)  integral  for’nula  'or  t“*e  electric  field  comonents  is  used  (’’ef,  1, 
'pMenpiv  f,  'o.  n~3ai)  where  ai^ain  contribtitions  fro"  aperture  ed^es  are 
•'O"  1  end  volume  current  and  charge  densities  are  assumed  to  be  zero. 

T*"'  '"ode  ir  assumed  to  be  propa'',atin''.  down  the  coaxial  line#  whose  curvature 
nepl^cted,  i.e.,  the  propagation  down  the  line  is  assumed  equivalent  to 
fiat  a1'  r  -  a  straight  coaxial  cable  ^lef*  1#  Appendix  IV,  Pa<:es  IV-6  and 
^■-7\ 

fn  this  sectioa,  Che  field  cononencs  will  be  formulated  if  the  cable 
•'"re  in  freo-spnee.  Tron  (T-3-<a,  b)  in  Appendix  F  of  "lef.  1,  v’C  have 

fi.r  t‘'<'  electric  and  naqnrt  ic  fields  T  and  at  the  point  ^  (characterized 
Hv  cvli^dricnl  coordinates  (0,  r.))  due  to  the  fields  on  a  cable  slot  located 

at  the  point  r',  (chnneterized  hy  cylindrical  coordinates  ((i,  b'),  where  -jL 

’■  s  the  cahic  radius  and  b*  is  the  vorticj!  CDor''in«nte  of  the  center  of  a 
Sint''; 


-  ^  h'-'o'r  x  ’•  £'i"> 
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M(o.  %  z)  =  jl  .iS"  *  ia"))C  -  (f  ^ 

AP" 


fG  -  rjV  :;(£■■) )'.'T, 


r;.l-b) 


*  1 

u^erc  t:ne  'Vnendencp  e  ''*  *  is  assuned  and  whf^re 
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Slot  area 
Creen*s  function 


Jk" 


!  =  aJ^  -  lAp  cos  G'  +  (z  -  b')^ 


O' 
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^ut'/r.rd  nornal  unit  v'jctor  at  slot 


-  A  U  AC  A  ^ 

^  +  X  ^  *  5  iil) 


basis  vectors 
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Per-’itt’vitv  oC  free  spree  =  i  10  1  f3rrrts/-olor 


=  Majtietic  neneabilit"  of  free  space  =  a"  henrjer./neter 


For  the  ■^'F.'I  node 


S(i')  =  t"  0..?-a) 

i(i')  =  3”"^x’> 

■  here  3"  Lr  the  enit  basis  vector  in  the  ilirection  in  the  ir.'i  >' e-priee'l 
coordinate  systen  In  which  the  slot  position  is  represented  (See  'ef.  "oction 
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2.1  Field  Components  from  a  Single  Slot 

Fron  (2.7-a,  substicuted  into  C2,3'-n,  h)  ani  with  the  aid  of  (2.3) 
we  obt>i2n  the  rectangular  field  components  at  a  point  r  •  (o,  t,  z)  from  a 
single  slot  at  the  point  *  (At  t  b*)  as  follows; 
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an.l  where  and  are  the  length  and  width  of  the  slot  respectively.  Ter 
the  p'.irpose  of  facilitating  integration  over  the  cable  it  is  convenient  to 
secarate  out  the  ii'-dependent  factors  in  each  terra  of  (2.d-a,  h,  cl  and 
'  2.h-a,  h,  c).  To  this  end,  we  note  that  p'  •  P  •*.  Q'  and  express  those  onuat.ions 
in  the  forms 
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T^o  for-'.T 

f2.10-a,  b,  cl  and  f2.1Ua,  b, 

,  c)  arc  part icularl '■ 

usefjjl 

for  nurnoTi^s  integratin'?  over  ill  t?'c  «lots  on  the  cabl*?,  as  vail  hecone 
aroiront  in  Section  2.2.  Tc  stu'l"  the  fields  fro-n  a  sinrle  slot,  the  lorns 
''J.S-a,  S,  r)  nni  b,  r)  are  less  cu>rberso*iie  and  tliereforc  easier 

to  'jse.  '’everCina  to  those  Cor'^s,  we  can  obtain  the  radial  (o-directed) 
and  i:*in’ithal  f '^-d irected )  components  ‘’ror  (2,Z~<i,  b)  and  b),  Sotc 

‘'irsr  t!'at  for  an/  vector  X  -  '^  +  ^  ^  where  ^  and  ^  are 

jr.lt  basis  vectors  in  the  c  and  >  directions  respectl'elv, 


in 
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It  was  shown  in  Ref.  1  (Paje  4.36,  Tqs.  (4.51^4.63)'  that 

•[(+)  ^  j(-)  ^  ^(0)  ^  j  i 

n  n  n  n 

''sins  (2.13)  in  (2.15-a,  b,  c)  and  (2.;6-a,  b,  c),  wc  ootain  the  aoi'roxiiate 
field  cotiponents 
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ipproxinatG  radial  and  aziciuthal  con’)onents  of  t!'.e  ^iolds  can  be 
obCfiineH  fron  {2,19-a,  b,  c)  nnd  ^2,20-a,  b,  '.’ith  the  aid  of  '2.17-a,  b,'. 
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3«  Calculations  of  Received  Sirnals 

rnlcuMtion  of  the  received  signal  consists  of  the  follo.tnn  steps: 
“'er'  1:  ''alculation  of  the  fields  from  the  cable  at  an  arbitrary  point 
'  as  :{  the  cable  were  in  infinite  free  space 
“to:  2:  Deterninntion  of  the  plane  wave  spectral  representation  of 
the  fields  calculated  in  Step  1 

Sr  O')  1:  evaluation  of  the  effects  of  the  presence  of  the  ground  on 
spectral  fields  decernired  in  Step  C 

■'tep  W  C-iordinate  transfornations  on  the  superposition  of  the  spectral 
calculated  in  Steps  2  and  3  in  or'^-^r  to  prepare  for  the  use  of  that 
■?  ;  ner’'  ''u t  ion  as  the  input  to  the  scattering  process.  The  transfornation 
ic  froo  the  eart.i  coordinate  systen  (x,  y,  z  or  P»  n,  z)  to  a  system  whose 
ori::.’n  at  t'^e  center  of  the  scatterer  and  whose  z-axis  is  parallel  to 

iiroction  of  the  propagation  vector  of  the  wave  incident  on  the  scatterer 
‘rep  "a’culation  of  the  fields  scattered  toward  an  arbitrary  point, 
usin''  *')o  ''ar^'cr  scattering  program 

'step  >:  "valuation  of  the  ground-reflected  scattered  field  components 
nt  un  arh’tr.irv  point,  again  using  the  Parber  scattering  nrogran 

7;  coordinate  transformation  of  the  scattered  field  component 
''lirert  nliis  tround-ref Iccted )  from  scatterer-centered  coordinates  to  ear^’h 
coord  1  r*atcs 

't^'r  :  Integration  ofthc  superposition  of  field  components  calculated 

n  'fcTC  2,  ■,  6  and  7  at  the  antenna  over  the  nlano-wava  spectral  space 

to  pr*  ‘•ICO  'c  total  field  at  the  antenna 

■"he  bac  ''round  analysis  and  discussion  pertaining  to  Steps  3,  S, 

1 ,  7  1  r  ! 


'f  coverod  in  Reference  1,  Sections  5,  6,  7,  3,  and  ID 


3.1  Plane-wave  Spectral  !^eoresentation  of  Fields 


The  representation  of  a  p.eneral  electric  or  nagnetic  field  vector  as 
a  superposition  of  plane  wave  fields,  known  as  the  "plane-wave  spectrun" 
of  the  field,  was  discussed  in  considerable  detail  in  !?ef.  1,  Sections  3 
and  4  and  Appendices  II  and  HI.  Referring  to  Appendix  II  of  Ref.  1  we  repeat 
here  the  Fourier  transform  pair  given  by  Fqs.  (11.31  and  (II. 4),  with  slight 
nodi f ications. 

First,  from  (II. 3)  the  field  vector  at  the  observation  point 

^  *  (o,  o,  7.)  is  expressed  in  terms  of  its  "plane  wave  spectrun’  by 


vfj)  =  V(£.  7) 


,2, 

li  3.  e  <e  V 


"-space  in 
2  2 

which  3.  *  • 
h  2 


2  «=  x^  +  yX  ■  horizontal  part  of  j 
2  =  unit  propaeatlon  vector  -  3^  ^  +  6_. 

'”v  i  '  horizontal  part  of 

h  *  *  V ^  *  tagniturie  of  horizontal  part  of  ^ 

^1  =  f\Jl  -  =  absolute  value  of  vertical  conporent  of  ^ 

=  tan  '  =  direction  annle  of  horizontal  component  of 

:  ’•3.  sin 


.  i"-  « 


d-3  =  (i:’  dS 

*  X  y 

V.(3,)  «  Portion  of  plane-wave  soectrun  of  V''r)  corresnondin''.  to 

^  "  +i!<|Sjz 

e  ‘  (upv/ard  propagating  plane  wave) 
a  Portion  of  plane-wave  spectrum  of  correspondinr,  to 


-jklEj 


(downward  propagating  plane  wave) 


Inversion  of  (3«2)  yields 


1#^  rj« 


The  work  that  was  done  on  this  project  to  calculate  olane-wave  spectral 
fields  5’(^)  fron  the  cable  fields  z)  and  — -  throu-’h  (3,2)  (an 

anpronch  that  was  eventually  abandoned:  see  Section  i)  is  rietaile''  In  Appemlix 
II.  Also  described  in  Appendix  Ti  is  the  work  done  on  an  alternative  to 
solution  of  (3.2)  used  to  obtain  an  aporoxination  to  the  pinne-wove  spectrun. 


J  J 


3>2  Ground  .^efloctions 


Reflection  of  a  plane  wave  with  wave  propagation  vector  propa^.atir'', 
in  seni-inCinite  free-space  bounc*od  oy  the  surface  of  the  earth  v/r.s  discussed 
in  Ref.  I,  Section  5-  The  results  are  ’iven  in  Tqs.  (5.13-n,  b).  “ith  srall 
chnnoes  in  notation,  the  results  for  the  electric  field  (na^^netic  field  result 
will  not  be  needed  here)  can  be  written  in  the  forn 

vher<? 


3.7 


3.3  Tields  Incident  on  tho  Scatterer -Coordinate  Transforr^ations 

coordinate  transformations  required  to  adaot  the  field  components 
to  the  ^.arher  scattering  program  are  treated  in  Section  6  of  heC.  1.  Figure 
on  Pope  3-2  of  ?cf.  I  is  repeated  here  (nunbered  Figure  3.3.1  in  the 
oresent  re^'ort)  to  illustrate  the  **leo  frame**  and  ’*body  frame"  coordinates. 

As  is  evident  fror  Figure  3.3.1  the  electric  field  vector  of  the  incident 
■nvo  :  in  the  (x^  y^)  plane  of  the  "*^b  frame.*'  'fhe  direction  of  the 
ohis  o*  t^o  "hoov  frame"  (which  isparallel  to  the  lon^  dimension  of  the  snher- 
oi'^nl  s'-ntterer)  is  defined  bv  spherical  angles  fO  ,  t  )  referenced  to  the 

p  p 

lib  •rn*"''.  Since  the  electric  field  h.ns  been  referenced  to  tbs  "ground  Crar’.e" 
ccori ' nates  y,  z),  it  v/as  necessnrv  to  develop  transformations  between 
the  groiio'^  frame  nnd  tie  lab  fra*‘e  and  also  bcfvoen  the  lab  frame  and  body 
,'r*-e.  These  tranTfornations  ••••“re  developed  in  Section  6  of  Ref.  1. 

■"’.e  transfornat ion  matrices  hot'^esn  body  and  lab  frame,  body  and  ground 
frame,  nrd  l.nh  and  ground  frame  (in  both  directions  for  each  case)  are  given 
bv  (A.3-r,  n',  (^.13-1,  b)*  and  (>.17-.n,  h)*  of  Ref.  I  respectively.  These 
^re  u'^e'  to  eve^op  the  expressions  for  the  electric  field  components  in 
the  1  ih  frn'-e  in  terms  of  those  in  the  '•round  frame.  "These  expressions  are 
MV?''  *•  V  ‘'•).2T-a,  0,  c)*  in  Rc^.  I  and  apply  to  plane  wave  spectre!  components 
<*f  \  v**  iel  '.  -rltter  in  matrix  form,  the  expressions  are 


1^.  1  .  I?.  1  «  Spectral  incii' 


field  vector  in  ’round  frnne  (3.'»-a) 


[3.  ’ 

*■  iz- 

I  Spectral  incident  field  vector  in 

1 


lab  frane  (3.^-b!'' 


L~.  i 

f'ote  that  r.  »  0,  as 


"cos  3g  cos  *2  cos  Oj 
-  sin  bp  sin  bj) 


from  Ti’ure  3.3. l-b' 


(cos  Sg  sin  bg  cos  bj  -  sin  3^ 

+  cos  bj  sin  <!j) 


I 

-  (cos  Og  cos  o,  sir  Oj  -  (cos  9g  sin  ij 
+  sin  bg  cos  bj)  -  cos  be,  cos  bj 

i 


sin  Oj  sin  0,  sin  bj 


I 

'_rin  cos  b^ 


sin  Gj  sin  bg  cos  3  , 


(3.W)' 

'ely  of 


vhere  anr  bg  arc  the  snltcrical  polar  and  azimuthal  angles  respectively  of 
the  wave  vector  <5  in  the  jround  frame  and  the  an’le  bj  is  defined  in  T.qs.  (f..''3-n, 
....  1  ')  of  bef.  I  ^or  the  tvo  scattering  processes  treated,  namely  scattering 
(iirectly  into  tl.e  antenna  ^Subscript  A)  and  scattering  toward  the  ground- 
ref'cction  point  f'libscript  b); 

cos  ^3  cos^bp  -  b^,.)  -  (Zg  -  Zj.)  Bin  Og 


point  f Subscript  b); 

Og,  cos  3p 


%  cos,bp  -  b^^)  -  (Zg  -  sin  Og 


I  (  beg  cos 


3,  COsf.g  -  hgg)  -  (Zg 


Zg)  sin  0,]^  +  [ogg  Sin(- 
(3.5-1) 


Sin('',  -  Sgg) 


vVh 


•'  '-•-V.-J 


V  <•  <*.  d.  d 


‘)etwpen  scatterer  center  and  antenna  and  j»round  reflection  point  resnect ivel y , 
and  f^qs.  0.5-e,  D*)  are  the  coordinates  of  the  Rround  reflection 


“he  required  Inputs  to  the  *^arber  scattering  program  are  the  incident 
':rld  o'^nlitude  ^Fq.  f^.33)  of  Fef.  I' 


Jr*,,  •  'Vi'  •  f*i,i 


inH  t''e  pr  1  .ir izat  1  on  an^le  of  the  incident  field  (nq.  n*  Ref  1  and 

'■cure  '.3.1  (h)  of  the  present  report) 


fhe  quintities  in  ''l.h-a,  b)  are  obtained  from  Rqs.  0.4i  with  the  ai'^ 

,,t  '3.R-a . c3',  (3. '-a,  ....  d)  and  (3.3-a . f)’. 

T^e  equations  '  fS-a,  h)  were  implesiented  on  the  conputer  as  a  part 
if  the  nrjeinal  project,  ’’'his  was  done  through  a  subroutine  called  "Rib." 

Ml  of  i-ie  Tortran  statements  in  Subroutine  RIS  and  the  calls  to  31?  from 
the  main  program  were  checked  as  a  part  of  the  current  project  and  it  was 
confirmee  tmt  a' 1  of  the  pro'jramming  -.as  f.mithful  to  the  analytical  results  summer 
izer  .roove.  3*,^  analysis  itself  was  also  checked  and  found  to  be  correct 
-ithi",  ‘he  ass'imptions  made  in  its  development. 


3<4  The  Scattering.  Process 

For  the  original  project,  the  scattering  process  was  modellefl  through 
the  scattering  program  developed  by  Peter  Barber.  This  program,  our  study 
of  alternative  possibilities  and  the  reasons  why  the  Barber  program  v.as  chosen 
as  the  means  of  modelling  scattering  from  a  human  frame  target,  are  discussed 
In  detail  in  Section  7  of  Reference  1,  The  theory  behind  the  program  is 
discussed  in  References  41,  45,  47,  50,  52,  57,  58,  60,  61. 

Further  investigation  of  alternatives  was  undertaken  at  the  beginning 
of  the  present  project.  The  possibility  of  changing  to  c  simple  short  Hioole 
scatteror  was  considered,  as  was  the  possibility  of  a  simplified  model  in 
which  the  scatterer  was  treated  as  a  small  perturbation  in  constitutive 
parameters  witun  the  volume  that  it  occupies.  The  snail  dipole  was  dismissed 
on  the  basis  that  the  human  frame  dimension  is  comparable  to  a  wavelength 
at  the  frequencies  of  interest  and  the  short  dipole  Is  an  extreme  oversimpli¬ 
fication.  The  second  idea  seemed  to  show  some  promise  In  reducing  computer 
time  (one  of  the  principal  weaknesses  of  the  Barber  program)  but  again  it 
is  not  as  accurate  as  the  Barber  program  and  the  limited  resorrees  at  our 
disposal  did  not  allow  the  extensive  software  changes  that  would  have  been 
necessary  to  implement  the  idea  for  this  application. 

for  the  above  reasons  it  was  finally  decided  to  continue  to  use 
the  Barber  program.  .\'',ain,  as  in  the  original  project,  it  was  implemented 
as  a  subroutine  called  by  the  main  program.  Some  snail  software  were  mode 
in  the  calling  format  within  the  main  program,  but  not  in  the  3nrbcr  subroutine 
Itself. 

The  Barber  nrogran  nodeis  the  sertterer  as  an  elcctroragneticall y  uniform 
spheroid  of  length  and  radius  where  and  2?^  arc  assigned  values 


corresponrti HP  to  the  height  and  :}irth  of  a  huinan  body.  The  constitutive  parameters 
for  the  scatterer  are  chosen  at  values  representing  a  composite  of  the  widely 
variable  constitutive  parameters  of  the  various  kinds  of  tissue  found  in 
thn  body  (e.p..,  bone,  skin,  fatty  tissue,  muscle,  blood).  These  values, 
cnrefullv  chosen  by  Professor  Barber  and  his  colleagues  in  connection  \/ith 
their  studies  of  absorption  of  RF  radiation  by  the  human  body,  ’vere  transmitted 
to  the  Principal  Investigator  on  our  project.  Those  were  the  values  used 
in  our  iTplenentation  of  BarberVs  program. 

The  Input  to  the  scatterer  in  this  proprar  is  assumed  to  be  a  linearly 
polari^e'^*  plane  i/ave.  That  is  one  of  the  major  reasons  for  our  use  of  the 
pi wave  spectral  representation  of  the  fields  Incident  on  the  scatterer 
^^ection  3.1  of  the  present  report).  The  approximation  described  in  “Section 
4,  vhich  wa«  eventually  invoked  in  our  final  computations  was  a  response 
to  the  difficulties  encountered  due  to  the  tise  of  the  plane  wave  spectrum 
concept.  Mtbough  the  scheme  described  in  Section  A  is  not  ideal,  it  is 
a  means  of  approximation  of  the  input  to  the  scatterer  that  captures  the 
most  imnortant  features  of  plane  wave  fields.  This  is  further  discusser 


3.S  Fields  Rmerains  from  the  Scattering 
Process-Coordinate  Transformation 


The  coordinate  transformations  required  for  the  scattered  fields  were 
described  in  Ref.  1,  Section  3.  The  geometry  is  illustrated  in  ^inure  ^ 
of  Ref.  1,  reproduced  in  the  present  report  as  Figure  3.5.1. 


Figure  3.5.1 
Scattering  Ceonetrv 
for  scattering  toward  antenna 

for  scattering  toward  ground-reflection  point 
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reflect  ion 


is  evident  from  the  figure,  the  scattered  field  is  neasured  alonr  the 
(X|  -  Zj^)  plane  in  the  lab  franc.  The  wave  vector  ‘^or  the  scattered  field 
■?  called  Tf.e  snnerlcal  polar  angle  of  in  the  la'  trice  is  denolr  ' 

cj  .  The  ontp'its  of  the  Barber  scattering  program  are: 


3 . 1 S 


M)  the  "vertically  polarized*'  scattered  field  component  i.e,. 


the  component  nornal  to  the  (xj  "  )  plane 


(2)  the  "horizontally  polarized"  scattered  field  component  i.e. 


the  component  parallel  to  the  (x^  -  z^)  plane. 


The  task  to  be  performed  is  to  obtain  the  scattered  field  components 
in  the  ’round  frame  from  the  components  in  the  lab  frame.  The  transformation, 
dcvploped  in  Section  8  of  Ref.  I,  is  given  by  Eqs.  (8.4-a,  b,  c)  of  Pcf.  1, 
repeated  below  in  matrix  form  with  slightly  modified  notation: 


where 
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Ground-frame  components  of  scattered  field  (3.7-a)' 
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i.ab-frane  components  of  scattered  field  (3.7-li)' 
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The  expressions  (3.7)  and  (3.7-3,  b,  c)'  are  evaluated  on  the  comnuter 


for  each  of  the  scattering  processes  "A”  and  "G,"  A  being  that  applicable 
to  scattering  directly  toward  the  antenna  and  G  being  that  directed  nt  the 
ground-reflection  point. 

The  relationship  between  the  lab-frame  components  of  F.g  and  the  horizontally 
and  vertically  polarized  components  and  is  given  in  (3.3-a,  b,  cl 
of  Ref.  1,  cast  in  matrix  form  as  follows: 


f3.E) 

where  is  given  by  (3.7-b)', 

where 


'•SH 
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SL 
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(3.B-n)’ 


and  where  cos  Ggj^  and  sin  0gj  are  given  through  (fi.5-a,  b)  and  ('.i— n,  b,  c 
of  Ref.  1. 

Based  on  (3.7),  (3.7-a,  b,  c)',  (3.8)  and  (3.8-b,  b)  of  the  'iresent 
report  and  F.qs.  (8.5-a,  b)  and  (8.6-a,  b,  c)  of  Ref.  1,  the  final  cor.rutr.tio" 
of  the  scattered  field  components  in  the  ground-frarre  in  ter-s  of  an'' 

can  be  achieved.  That  result  is  given  in  fb.|l_a,  p,  c)  of  Ref.  1,  It  is 


reneated  belo.'  in  matrix  form; 


fEi  '  and  [fi  1  are  given  by  (3.7-a)*  and  (3.3-a>'  respectively  and 
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where 


(5.G-a)' 


cos  and  sin  4j  for  both  A  and  G  processes  are  given  by  d3,5-a,  d) 
wit*'  the  aid  of  (3.3-a,  f)* 


cos  0, 
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P.ef.  1) 

,(U) 

’si. 

’’su 
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l*T  j|l 
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where 
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‘'su 


rectang'jlnr  coordinates  of  point  I'  (3.9-d)' 

I  2  2  2 

^\r’'u  "  ’'s^  +  (y;-  -  y^) 

between  sentterer  center  and  point  U  (3.9- e)' 
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3.18 


z,  j,  -  sin  Gg  cos  «ig  (x,|  -  Xg)  +  sin  Gg  sin  <3^  (y|.  -  Vg) 

+  cos  Gg  (z,j  -  Zg) 

and  where  I.'  is  either  the  antenna  (Point  A)  or  the  nrounc!  reflection  noint 
(Point  C). 

Tn  the  original  project,  the  equations  above  were  inplenenled  in  n  sub¬ 
routine  called  ''50S."  The  analysis  on  which  BOS  is  based,  the  '^ortrnn  in 
the  subroutine  itself  and  the  calls  to  it  fron  the  main  pro*’, ran  \.'cre  all 
examined  thoroughly  as  a  part  of  the  current  project.  The  analysis  was  found 
to  be  correct  within  the  assumptions  made  and  the  programning  was  found  to 
he  consistent  with  the  results  of  the  analysis,  so  no  ch.anges  were  made  in 
the  software  relating  to  the  coordinate  transformations  at  the  output  of 
the  scattering  process. 


I.", 


3.10 


3.6  rffects  of  Groun'*.  Pe^'Iections  on  the  Scattered  Fields 


Pie  effects  of  ground  reflections  on  the  scattered  fields  was  covered 
in  'action  °  of  Ref.  1.  It  is  an  adaptation  of  the  theory  discussed 

in  Section  3.2  of  the  present  report  to  the  special  case  of  a  plane  wave 
prQna'’atinp,  in  the  direction  from  the  scatterer  center  toward  the  ground 
reflection  point.  Although  the  scattered  wave  is  actually  a  spherical  wave, 
the  groun'i-ref lect ion  process  is  rodelled  as  if  a  plane  wave  were  incident 
on  the  ground  surface  in  a  region  containing  the  ground  reflection  point 
'i.e.,  Che  point  at  which  the  law  of  reflection  Is  obeyed  for  the  path  frost 
scatterer  center  to  ground  reflection  point  to  antenna). 

The  wave  vectors  directed  from  the  scatterer  center  S  toward  the  ground- 
reflection  point  G,  denoted  by  is  given  by  Eqs.  (9.12-b,  c,  e)  in  Ref.  1, 

The  wave  vector  directed  frora  the  point  G  toward  the  antenna  A,  denoted  by 
is  given  by-"qs.  (O.U-f,  g,  1)  in  Ref.  1.  The  quantities  Ygj,  and 
for  these  wave  directions  are  given  by  {0.12-j,  k)  in  Ref.  1. 

The  following  equations,  taken  from  Ref.  1  (Eqs.  9.12-a,  ....  k),  sunmarize 
the  results  of  the  analysis  pertaining  to  ground  reflection  of  the  scattered 
fields: 
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Coordinate  of  ground 
reflection  point  relative 
to  scatterer  center  in 
terms  of  antenna  and 
scatterer  coordinates 
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distance  between  scatterer  center  and  ground  reflection  point 
in  terms  of  antenna  and  scatterer  coordinates  f3.10-d> 
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components  of  wave 
vector  for  scattered  wnve 


incident  on  ground  reflection  (3.ir>-h) 
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Parameters  for  ground-reflection  point  to’  antenna  nroonn.at icn 
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A  ,  ,  reflection  point  in  terms  .. 
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1  of  antenna  and  scatterer 
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Components  of  wave  vector 
for  scattered  wave 


reflected  from  ground  and 


>  propagating  toward  antenna 
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Civen  Cqs.  ^3.10-a,  ....  j)  and  (3.11-a,  j),  we  can  apply  the  theory 

pf  ground  reflection  of  plane  waves  summarized  in  Section  3.2.  "e  invoke 
'■.q.  0.3)  in  the  form  given  hy  f9.13)  in  Ref.  1.  The  final  expression  for 
the  ground-reflected  scattered  field  vector  at  the  antenna  is: 

(3.12; 


,7(GA), 

'  'S  ^  (  Sx 


(3.12-aV 
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as  given  by  (3.9)  for  the  (3.12-b)’ 
"G"  SC'  ►ering  process 


Eik  o'=<hn'  Ejk 


f3.12-c)' 


3.7  Total  Field  at  the  Antenna 


intesration  of  tie  components  of  the  spectral  field  over  spectral 
space  (i.e.,  the  space  of  the  vave  vector  2)  to  produce  the  components  of 
the  total  (direct  from  cable  plus  p,round-ref lected  plus  scattered  directly 
toward  antenna  plus  ground-reflected  scattered)  field  was  described  in  detail 
Ir,  Reference  1,  Section  10.  "Hiis  operation,  which  is  actually  a  two-dimensional 
inverse  Fourier  transformation,  is  given  genetically  by  Rq.  (10.1)  in  Section 
1C  of  Reference  1,  as  follows: 


(8-spacc) 

•/here*?(;^)  is  a  component  (e.g.,  x,  y  or  z)  of  the  total  spectral  field 
and  F,(]^)  is  the  same  component  of  the  actual  total  field  as  a  function  of 
the  position  coordinates  of  the  observation  point,  is  the  free-space 
wave  number  and  d  denotes  dfi^  dS^.  Using  the  cylindrical  coordinate  repre¬ 
sentation  of  ^  (i.e.,  g  •  3^(5  cos  6g  +  J  sin  *g)  ±  |3^|^)  and  noting  that 
2 

•\1  -  .3^  ,  the  Integral  (3.13)  can  be  converted  to  the  form 


"(0,  di,  z) 


Jk^[p6j^cos(6g  -  0)  z^l  -  3 


(3.16) 


vhere  it  is  noted  that  only  the  "upgoing"  spectral  wave  (i.e.,  ^Sj,l) 

enters  Che  computation,  because  both  direct  (from  cablel  and  ground-reflected 
spectral  waves  incident  on  the  scatterer  are  iipgoing. 

Tbe  procedure  for  integrating  (3.14)  in  the  original  project  was  to 
first  evaluate  the  integral  numerically  and  then  evaluate  the  dg  integral 
bp  stationary  phase  methods.  It  was  recognized  that  a  compromise  in  accuracy 
was  incurred  in  certain  parameter  regimes  through  the  use  of  stationary  phase. 
However,  the  expenditure  of  computer  time  in  performing  a  more  rigorous 


nunerical  integration  on  is  prohibitive.  Consequently  it  was  deciiled 
to  Continue  to  use  stationary  phase  for  the  Ag-integratlon  in  the  present 
project. 

There  was  some  improvement  in  the  technique  for  carrying  out  this  double 
integration,  however.  A  new  numerical  integration  procedure  was  adopted 
for  the  3^  integral  that  is  very  fast  and  requires  only  12  points  to  attain 
great  accuracy  over  the  range  from  3^^  •  0  to  3^^  •  1.  That  technique  was 
used  in  the  computations  that  were  done  using  the  plane-wave  spectrum  approach. 
In  Section  4,  a  new  ap.oroach  not  requiring  this  integration  Is  described, 
and  the  latter  approach  was  used  to  obtain  our  final  numerical  results. 

Tence  the  issue  of  how  the  integration  is  performed  becomes  academic  with 
respect  to  the  actual  results  presented  in  Section  5. 
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4.  Approximation  of  Field  at  Scatterer  as  a  Single  Plane  ''ave 

In  Appendix  II  the  difficulties  experienced  in  trying  to  inpleraent  Che 
plme-'^ave  spectrum  aporoach  are  cescribed  in  some  detail.  In  response  to 
these  difficulties,  a  new  scheme  was  devised,  in  which  tne  plane-wave  spectrum 
approach  was  abandoned  in  favor  of  an  approach  wherein  the  field  incident 
on  the  scatterer  is  approxirruated  as  a  single  plane  wave  (i.e.,  not  as  a  super¬ 
position  ot  plane  waves)  Obviously  this  constitutes  a  sacrifice  in  realism, 
bee  luse  the  field  Incident  on  the  scatterer  from  the  cable  is  not  a  plane 
wavs.  However,  the  scattering  is  deterntned  by  the  incident  field  within 
the  volume  occupied  by  the  scatterer.  The  incident  field  pattern  outside 
the  scatterer  volume  does  not  affect  the  scattered  wave.  Thus  if  the  varia¬ 
tion  of  the  incident  field  resembles  that  of  a  plane  wave  coming  from  a  particular 
direction  within  the  scatterer  volume,  regardless  of  its  behavior  outside 
th.et  volume,  then  it  can  be  so  approximated  for  purposes  of  determining  the 
scattered  wave  fields. 

Various  ways  of  generating  a  plane  wave  approximation  at  the  scatterer 
were  contemplated.  The  major  problem  was  to  determine  a  proper  direction 
for  the  incoming  "approximate  plane  wave."  It  was  decided  that  the  optimal 
direction  to  assume  was  that  in  which  the  power  flows,  i.e,,  the  direction 
of  the  Poynting  vector  for  the  total  field  from  the  cable,  i.e.,  the  fields 
whi se  rectangular  components  are  given  by  (2.19-a,  b,  c)  and  (2.20-a,  b,  c). 

The  procedure  was  to 


(1)  Calculate  the  Poynting  vector  for  the  total  E  and  H  fields  from 
the  cable  at  an  arbitrary  observation  point  (p,  p,  z).  This  complication 


(4.1) 


or  in  component  form 


-  iRe  (Ej,H*  -  E^l!*)  (4.1-a)' 

where  E^,  E^,  F.^,  H^,  and  arc  given  by  (2.19-a),  (2.19-b),  (2.19-c), 
(2.20-a),  (2.20-b)  and  (2.20-c)  respectively. 

(2)  Calculate  the  amplitude  of  the  Poynting  vector 

Pfl  •  SJP^  +  Py  +  Pj  (4-2) 

where  P^^,  P^  are  obtained  from  (4.1-a,  b,  c)*. 

(3)  issune  the  normalized  propagation  vector  ^  for  the  assumed  plane- 
wave  incident  on  the  scatterer  to  be  in  the  direction  of  the  Poynting  vector 

at  the  center  of  the  scatterer.  Since  by  definition,  must  h.>'v  unit  magnitude, 

this  implies  that  the  components  of  ^  are  given  by 

P 

3^  »  at  scatterer  center  (4.3-a) 


P 

S  •  at  scatterer  center  (4.3-b) 

P 

3  ■  at  scatterer  center  (4.3-c) 

o 

The  calculations  of  the  components  of  ^  are  made  from  (4,3-a,  b,  c) 

with  the  components  of  P  calculated  from  (4.1-a,  b,  c)'  and  P  from  ^4.2). 

o 

(4)  The  values  of  3  ,  £  ,  3  obtained  from  (4.3-a,  b,  c)  are  used  as 
X  y  2 

inputs  to  the  Sarber  scattering  program,  in  order  to  define  the  direction 
o;  the  plane-wave  assumed  to  be  Incident  on  the  scatterer. 


It  should  he  noted  that  the  approxination  used  here  is  based  on  the 
ieterminatlon  of  the  direction  of  power  flow  only  at  the  exact  center  of 


the  scatterer.  The  actual  direction  of  the  Poyntlng  vector  will  be  different 
It  different  points  on  the  scatterer,  but  we  are  in  effect  assuming  that 
the  Incident  wave  can  he  approximated  by  a  linearly  polarized  plane  electro¬ 
magnetic  wave. 

The  properties  of  the  assured  plane  wave  are  as  follows: 

(1)  the  amplitude  of  a  component  of  the  electric  field  of  the  incident 
plane  wave  is  Che  amplitude  of  that  component  of  the  actual  electric 
field  at  the  center  of  the  scatterer 
;2)  the  incident  plane  wave  is  propagating  in  the  direction  of  the 
Poyating  vector  (i.e.,  direction  of  power  flow)  of  the  actual 
fields  at  the  center  of  the  scatterer 
f3)  the  spatial  direction  of  the  electric  field  of  the  Incident  plane 
wave  is  the  spatial  direction  of  the  actual  electric  field  at  the 
center  of  the  scatterer 

Assumption  (2)  guarantees  that  the  assumed  plane  wave  will  be  transverse, 
because  it  forces  the  propagation  direction  as  defined  by  6  to  be  .normal 
to  both  electric  and  magnetic  field  vectors  at  the  center  of  the  scatterer. 
Although  the  phase  fronts  of  the  field  may  have  some  significant  curvature 
within  Che  scatterer  volume  if  that  volume  is  comparable  to  a  cubic  wavelength, 
and  the  v.nrlatlon  of  phase  in  the  assumed  propagation  direction  ^  is  not 
ex.nctly  the  same  os  that  of  a  plane  wave,  the  approxiTct?  wave  field  is  that 
of  a  transverse  linearly  polarized  wave  in  the  region  near  the  center  of 
the  scatterer.  To  repeat  the  point  made  above  about  the  nature  of  the  assumed 
plane  wave,  if  the  propagation  direction  as  defined  by  ^  is  parallel  to 
the  Poynting  vector  at  the  scatterer  center,  then  ^  must  be  normal  to  both  the 
electric  and  magnetic  field  vectors  at  the  scatterer  center,  guaranteeing 
th-nt  the  .issumed  plane  wave  is  transverse  throughout  the  volume  of  the  scatterer. 


It  is  also  linearly  polarized  in  that  resion,  because  the  direction  of  the 
electric  field  of  the  incident  wave  will  be  assumed  to  be  that  at  the  scatterer 
center.  That  direction  is  fixed  as  that  o*^  the  actual  total  electric  field 
f roe:  the  cable  at  that  point. 

Software  was  written  to  implement  the  computations  of  P  ,  ^  ,  P  Cron 

r  .  X  y  ’  a 

(4.1-a,  b,  c),  (2.19-a,  b,  c)  and  (2.20-a,  b,  c),  the  subsequent  computation 

of  P  from  (4.2)  and  computation  of  3  ,  S  ,  3  from  (4.3-a,  b,  c'.  The  end 
o  x  y  z  V... 

products  of  these  computations  are  the  components  of  which  define  a  propa¬ 
gation  direction  of  the  assumed  plane  wave  incident  on  the  scatterer  an! 
which  are  inputted  into  the  Barber  scattering  program.  The  electric  field 
amplitude  and  direction  (in  a  plane  normal  to  g)  for  that  plane  wave,  assumed 
to  be  those  of  the  actual  computed  field  at  the  scatterer  center,  are  also 
delivered  to  the  Barber  scattering  program  as  the  amplitude  and  polar’ zation 
direction  of  the  incident  plane  wave.  The  phase  of  the  incident  plane  wave’s 
electric  field  at  the  scatterer  center  is  assumed  to  be  the  sne  as  that 
of  the  actual  field  at  the  scatterer  center  and  to  vary  along  the  propagation 
direction  (i.e.,  in  the  direction  parallel  to  as  docs  that  of  a  plane 
wave  Propagating  in  the  direction  of 

■dhis  plane  wave  approximation  is  assumed  to  hold  within  the  scatterer 
region  for  both  the  direct  field  from  the  cable  and  the  ground-reflected 
field,  ’^he  latter  is  also  approximated  ns  a  plane  wave  and  its  'lireci’on 
within  the  scnttcrer  volume  is  assumed  to  be  the  same  as  that  of  the  ’.ave 
directly  from  the  cable.  This  is  a  valid  approximation,  because  the  source 
of  the  field  is  (for  all  practical  purposes)  at  ground  level  and  hence  the 
ground  reflection  point  can  be  located  at  the  source  with  negligib’e  error. 

The  elevation  of  the  cable  slots  rolati''c  to  ground  level  Is  b'  .Cy'TO, 
less  than  1  centimeter  above  ground;  hence  for  the  scatterer  center  nt  its 
lowest  possible  point  within  the  parameter  regime  of  interest  in  f"!'  stu ■ 


(e.!'.,  the  case  of  Che  crawling  intruder,  where  the  height  of  Che  scatterer 
center  is  2n  c-i.),  the  ground  reflection  point  is  extrenely  close  to  the 


cable.  To  show  this  note  that  the  distance  between  the  cable  and  the  ground 
reflection  point  is  (See  Figure  4.1) 


ohere  d  is  the  radial  distance  between  the  cable  and  the  scatterer  center, 
bj  is  the  height  of  a  cable  slot  and  h,  is  the  height  of  the  scatterer  center. 
For  the  crawling  Intruder  at  the  furthest  radial  distance  in  our  computations 
(where  d  -  6.5  m,  h^  -  .0069  n,  h,  •  .25  si),  the  radial  distance  from  the 
rahle  to  the  ground-reflection  point  Is  only  about  22.5  cm  and  it  is  still 
srallor  for  larger  scatterer  heights  and  smaller  scatterer  distances.  To 
put  this  another  nay,  the  ground-reflected  wave  at  the  scatterer  center  appears 
to  originate  at  the  inage-point  shown  in  Figure  4.1,  which  isahout  .69  cm. 
below  ground  level.  The  tangents  of  the  elevation  angles  of  ^  and  (the 
amles  F  and  respectively  in  Figure  4.1),  in  the  "worst  case"  alluded 
to  above,  are 


.0669  rdn 


.0398  rdn 


2.115  degrees 

2.231  degrees 


f4.5-a) 


(4.5-h) 


■Fne  lif^'ermcc  in  angle  between  the  propagation  vectors  of  the  direct  and 
groiind-reriectod  ^assumed)  plane-waves  in  the  worst  case  is  only  .166  degrees, 
hrrcp  for  poactical  purposes  the  propagation  directions  for  these  two  waves 
can  he  ronsi^'ered  parallel. 

■’.ased  m  the  aliove  arguments,  the  ground-reflected  wave  impinging  on 
f-e  scatterer  is  modelled  as  a  plane  wave  propaeating  in  the  direction  of 
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the  Poynting  vector  for  the  total  field  from  the  cable  at  the  scatterer  center. 
The  amplitude  and  phase  of  the  components  of  that  wave  at  the  scatterer  center 
are  those  obtained  through  the  ground-reflection  transformation  epuations 
(3.3)  and  (3.3-a,  b,  ....  f,  g)',  where  the  parameters  (3|^,  dig)  used  in  (3.3-c, 
d,  e,  f)'  are  otained  from  (3_.  3„,  8.)  as  calculated  from  (4.3-a,  b,  c). 


f!°3ults  and  Conclusions 


'n  thi<;  section  the  nunerlcnl  results  are  prosentch  C^ubsection  5.1) 
in'^  then  discussed  (Subsection  5.2). 

5.1  N'unerlcsl  results 

b'unerical  computations  were  nade  of: 

M)  The  rectangular  components  of  the  (direct  and  ground-reflected) 
electric  fields  at  various  observation  points  near  the  cable  in 
the  absence  of  a  scattercr  (i.e.,  the  '’llluininating  fields"  in 
the  region  where  the  scatterer  would  be  located) 

(2)  The  rectangular  components  of  the  (direct  plus  ground-reflected) 
electric  field  at  the  antenna  In  the  absence  of  a  scatterer 

(3)  The  rectangular  components  of  the  contributions  to  the  field  at 
the  antenna  due  to  the  presence  of  the  scatterer  (wave  scattered 
directly  into  the  antenna  plds  ground-reflected  scattered  wave 
at  the  antenna) 

C41  The  rectangular  components  of  the  total  field  (direct-wave-pliis- 

ground  reflected  wave  plus  direct  scattered  wave  plus  ground-reflected 
scattered  wave)  at  the  antenna 

The  format  for  these  computations  was  essentially  the  same  as  that  in 
'Reference  I,  Section  H.  For  the  set  of  computations  (1),  the  amplitudes 
of  the  rectangular  component  of  the  electric  field  were  computed  for  eight 
values  of  0,  the  radial  coordinate  at  the  observation  point,  at  a  fixed  value 
of  the  azimuthal  coordinate  0  and  the  height  z.  Tor  one  particular  set  of 
results  presented,  z  is  set  at  0.5  meters  and  b  is  set  at  O’,  100°  or 


270°.  "or  mother  sot  of  results  presented,  <>  is  set  at  0'  and  a  is  set 
at  l.a,  2  anr  2.5  meters,  "or  each  fixed  pair  of  values  of  ■  and  z,  a 


5.2 


set  of  runs  was  made  for  P  >  17.5,  19.0,  20.5,  22.0,  23.5,  25.0,  26.5  and 
28.0.  In  each  case,  the  quantities  |K^| ,  IK^I ,  I  I  |E^|^  and 

IeJju  -  20  log^Q  lEj,  iE^I^  and  were  computed. 

The  conputations  (2)  were  also  of  the  absolute  amplitudes,  squared  ampli¬ 
tudes  and  amplitude?  in  decibels  of  the  rectangular  components  of  the  electric 
field,  but  at  a  single  point,  the  location  of  the  antenna.  The  assumed  antenna 
coordinates  are  o  =  0  and  z  •  1. 

The  forr^at  for  the  computations  (3)  and  (4)  was  exactly  the  sore  ns 
that  for  the  computations  (1).  Of  course,  the  coordinate  values  p,  "i  and 
z,  which  were  the  coordinates  of  Che  observation  point  in  computations  (1), 
are  now  the  coordinates  of  the  center  of  the  scatterer.  The  coordinates 
of  the  observation  point  in  conputations  (3)  and  (4)  are  those  assumed  for 
the  antenna, D  ■  0  and  z  ■  1. 

The  purpose  of  the  computations  (1)  is  to  determine  the  amplitudes  of  the  components 
of  the  fields  that  illuminate  the  scatterer  as  its  distance  from  the  cable 
varies.  The  cable  is  at  p  approximately  equal  to  24  meters;  hence;the 
radial  distance  of  the  observation  point  from  the  cable  is  varied  from  6.5 
meters  inside  (p  «  17,5)  to  4.0  meters  outside  (p  •  28). 

Computations  (2)  and  (3)  were  done  for  the  purpose  of  separatin"  out 
Che  contributions  to  the  fields  at  the  antenna  location  from  the  cable  (both 
direct  and  ground-reflected)  and  tne  contributions  from  the  scatterer,  Compu- 
tations(4)  are  those  of  the  end-product  of  this  study,  the  total  field  due 
to  cable,  ground  reflections  and  scatterer. 

Since  the  antenna  used  in  t'le  system  of  lircct  interest  is  vertical 
and  therefore  responds  only  to  the  vertical  electric  field  component,  plots 
are  presented  only  for  the  vertical  field  components  in  Computations  (3). 

Tabulated  results  ere  prpsented  for  .all  field  components. 


,  w.  . 


,  ■ 
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For  the  computations  (3)  and  (•.),  three  cases  are  considered,  as  in 
Reference  1.  The  first  case  is  the  "radial  walk,"  where  the  spheroid  that 
represents  the  human  frame  target  has  its  long  axis  in  the  vertical  direction 
as  it  noves  between  O  -  17.5  and  o  »  28.0.  The  second  case  is  the  "radial 
cravl,"  where  the  spheroid's  axis  is  in  the  horizontal  direction  and  pointed 
in  the  radial  direction  as  it  traverses  the  pcth  between  P  17.5  and  p  =  23. 0, 

The  third  case  is  tho  "radial  walk  on  stilts,"  where  the  upright  ""lan"  is 
elevated  hy  1  meter  relative  to  ground  level,  implying  the  use  of  1  mctcr-iii’h 
stilts  to  elude  detection. 

The  tables  below  (Tables  5.1,  5.2,  5.3)  constitute  an  account  of  the 
numerical  results  for  Computations  (1),  (3)  and  (4).  The  results  are  all 
presented  in  decibels  relative  to  a  reference  level  of  1  v/m.  In  each  column, 

Che  peak  value  (If  It  exists)  is  Indicated  by  surrounding  Che  number  by  a  rectangle.* 
The  results  for  Computation  (2),  the  fields  at  the  antenna  without  the  scatterer, 
are: 
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*■  In  some  cases,  there  is  no  clearly  definable  peak  value,  in  which  case  there 
'S  no  such  indication. 
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5.2.1  The  illuninatina  fields 


Consider  first  the  coraponents  of  the  "illuminating  electric  field  with 
no  ground  reflections,  as  given  in  Table  5.1  and  Figures  5.1-5.12.  These 
results  are  fron  conputations  of  the  amplitudes  of  the  quantities  in  Tqs. 
(2.19-a,  b,  c).  Since  these  results  involve  the  integrals  around  the  cable, 
ve  cannot  say  that  the  fields  originate  from  any  ore  particular  portion  of 
the  cable.  However,  intuitively  (aided  by  the  discussions  in  Appendix  I, 
expecially  Section  I-C)  it  is  evident  that  at  points  very  close  to  tie  cable, 
the  major  contribution  to  the  fields  should  come  fron  the  slots  closest  to 
the  observation  point,  i.e,,  those  slots  where di'  is  very  close  to  b,  'Isinj 

V'“'y  "  2 

A"  +  (z  -  b')  (where 

A  •  -  0)  the  d  istance  from  the  cable  slot  corresponding  to  ■  b,  should 

be  easy  to  establish.  First,  the  factors  (J!:  -  Ij)  that  appear  In  all 

the  electric  field  components  of  (2.19-a,  b,  c)  can  be  thought  of  ns  varying 


roughly  as  Ari  l  *  -■  .  Tn  E  and  E  ,  there  is  a  factor  proportional  to  A 

6^  V  My  *  ^ 

multiplying  (jk  Ij  -  I^)  in  one  of  the  terms.  Thus  the  factors  A(jk  -  Ij) 
in  the  cos  b  and  sin  0  terms  in  E  and  E  respectively  should  vary  roughly 

*  y 


-T\|*  ♦ 


The  sin  *  term  of  E  and  the  cos  ®  tern  of  E  should  vary 
“  ) 


roughly  as  Based  on  these  considerations*  ve  should  be  able  to  conclude 
from  (2.19-a,  b,  c)  that 


At  *  -  0®  or  *  -  18C“ 


r  ^r\li  +  — 
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.,  f)  tend 

to  corroborate 

(5.2-d) 


(‘5.2-e) 


f5.2-f) 


The  najor  trenHa  observed  in  Table  5.1,  Parts  (a),  (bi,  (c)  ere: 

('. )  The  peak  values  of  all  three  field  amplitudes  always  occur  at 
p  ■  23.5  m. 

(2)  At  P  «  23.5,  where  peal:  values  occur,  the  y-components  of  the  field 
amplitudes  at  P  •  00°,  ISO®  and  270°  are  nearly  the  same  as  the 
x-components  at  4  =  0°,  90°  and  180°  respectively,  l.e.,  the  x 

and  y  components  appear  to  interchange  their  roles  at  90°  intervals. 

(3)  At  other  values  of  p,  the  x  and  y  field  amplitudes  show  very  weak 
dependence  on  4. 

'4)  At  0  =  23.5,  peaks  of  the  x-field  amplitude  are  roughly  the  snne 
•■’s  those  of  the  x  and  y  components. 

''5)  At  other  values  of  0,  the  x-component  values  are  considerably  lower 
tiian  the  x  and  y  component  values. 


5ome  of  these  effects  are  discussed  more,  quantitatively  below.  Tron 
2-a,  ...,  f),  wr  conclude  that,  if  cable  attenuation  is  neglected 
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Tho  ratio  %  V.R  in  ~able  5.1  is  about  -5  dD  at  D  =  23.5,  as  connared 
with  -3.5  d3  as  in  the  tabulation  below.  The  computational  results  p.iven  in 
Table  5.1  for  (1\j  at  p  «  25.0,  (20.5  or  26.5)  and  17.5  are  (-3  to  -4  dF), 

(j  to  -2  dF)  and  (-1  to  -3  dB)  respectively.  The  corresponding  values  on 
tbe  tabulation  above  are  -.71  dB,  -1.36  d?.  and  -2.19  d3  respectively. 

Tne  corputed  values  of  the  ratio  °  *  23.5,  25.0,  (20. 5  or 

26.5)  and  17.5  respectively  are  (0  to  -1  dO),  (-5  to  -10  dB),  (-14  to  -17  dF) 
and  (-23  to  -25  dQ)  respectively.  The  corresponding  values  on  tbe  tabulation 
above  are  0  dB,  -6  d3,  -15.6  dP,  and  -21.6  dB  respectively. 

For  tho  most  part  (l.e.,  e.xcept  for  the  resultyj^  at  o  =  25.0),  these 
results  indicate  that  the  illumination  of  the  scatterer  situated  within  about 
6  -eters  radial  distance  from  the  cable  on  either  side  are  primarily  due  to 
the  slots  near  tbe  scatterer  location,  i.e.,  the  rough  estinates  culnlnating 
in  the  tabulation  above,  which  are  based  on  the  approximation  that  the  field 
comes  entirelv  from  the  cable  slot  at  o’  s  m,  come  remarkably  close  to  the  results 
piven  in  Table  5,1,  (a),  (.b),  (c)  which  take  account  of  contributions  from 
all  portions  of  the  cable. 

The  remaining  dependence  on  m  in  (2.19-a,  b,  c)  arises  through  the  cable 
attenuation  factor  which  is  present  in  C(b).  This  number  In  dB  is 

-o.oB  C'A.'.  Ttb  the  values  ft  «  .G02,  A  >  24,  the  ilecibel  difference  at  90° 
intervals  "iie  to  this  factor  is  about  .417  ^  or  .65  dB,  Thus  the  variation 
with  i  due  to  this  attenuation  factor  over  the  entire  cable  is  about  2,6  dB, 
correspondin"  to  a  130°  variation  of  about  1.3  dB.  That  helps  explain  the 
snr.'l  decrease  i between  1  and  2  dB)  in  field  amplitudes  between  0°  and  130° 
m  l  betweer  90"  and,  270°.  Otherwise  the  variations  between  the  |F^.!  and 
iE,.l  resu'ts  ior  P  -  (0°,  130°)  and  t'losf  for  o  =  (9C°,  270°)  appear  to  be 
d”.e  to  t  'lft  '"ictors  cos  6  and  sin  o  in  ^2.19— e)  and  (2,l9-h),  Because  or 
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these  factors,  the  computed  values  of  at  ■■)  •  X“,  lf?0“  and  270“  are 

nearly  indentical  with  those  of  at  0°,  90®  and  180®  respectively. 

From  Table  5.1  (c),  it  is  evident  that  there  is  no  significant  dependence 
of  |C^|  on  o  except  for  the  saall  (1-2  dP)  difference  between  values  of 
C  separated  by  180®  due  to  the  attenuation  factor  alluded  to  above.  Othen  ise, 
as  indicated  by  the  absence. of  ^-dependence  in  Fq.  {2.19-c),  the  amplitude 
has  nearly  perfect  cylindrical  synmetry. 

'■'e  now  consider  the  variation  of  the  illuninating  fields  with  the  height 
of  the  observation  point  a  at  rt  «  0®  (Table  5.1  (d),  (e),  (f)).  As  a  varies 
fron  1  to  2.5,  the  peaks  gradually  disappear  and  in  general  the  variation 
of  field  c.nplitude  becoites  Increasingly  flattened.  This  is  clearly  due  to 
the  fact  that  the  dependence  of  field  strength  on  horizontal  distnnee  fron 
the  cable  is  weakened  as  the  vertical  distance  Increases.  Tiie  mathematical 
denendence  of  field  strength  on  these  two  components  of  the  separation  distance 
is  contained  in  the  distance  which  is  essentially  the  square-rcot  of  the 
sun  of  squares  of  horizontal  and  vertical  distance  from  the  cable.  As  the 
vertical  distance  is  Increased,  the  contribution  fron  the  horizontal  distance 
is  increasingly  "swanped  out"  by  that  involving  the  vertical  distapce. 

Finally,  we  observe  by  comparing  Parts  (a)  and  (b)  of  Table  5.1 

with  ’’art  (c)  of  that  sane  table  that  there  is  appreciable  difference  between 
the  vertlc.nl  and  horizontal  field  amplitudes  in  the  regions  near  the  peak 
values  but  that  the  vertical  amplitudes  decay  more  rapidly  as  the  observation 
point  recedes  from  the  cable.  The  reason  is  easily  reternined  from  Fqs.  (2,19- 

3,  j,  c).  The  vertical  field  decays  roughly  es  and  the  horizontal  fields 

1  - 
as  -j  as  the  distance  from  the  cable  increases.  Thus,  while  their  values 

are  comparable  near  the  cab’e,  the  vertical  fields  decay  more  ranldly  with 

distance  from,  the  cable. 
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5.2.2  The  scattered  fieltls 

'.'e  nov;  turn  to  the  results  for  the  scattering  contrihution  to  the  field 
conponents  at  the  antenna  (Table  5.2  rnd  Figures  5.13-5.24). 

First,  've  consider  the  variation  «ith  o  for  fixed  i  and  z.  There  are 
throe  sources  of  this  variation,  as  follows: 

(1)  The  variation  of  the  fields  illuninating  the  scattercr  with  the 
distance  from  the  cable.  Tliis  was  discussed  in  Section  5.1  and 
is  surmarized  in  Egs.  (5.2-a,  ...,  f).  If  the  effect  of  ground 
reflections  on  the  illuninating  fields  is  included,  these  results 


will  be  changed  but  the  variation  with  distance  fron  the  cable 
will  not  be  radically  different. 

(2)  The  factor  — ^  on  the  wave  scattered  directly  into  the  antenna 

and  the  factor  /'  —  - r  on  the  ground-reflected  scattered  wave, 

*  ‘'ga' 

where  tg^  is  the  distance  between  scattercr  center  and  antenna 
location,  r^^  is  the  distance  between  ground-reflection  point  and 
antenna,  and  is  the  distance  between  scatterer  and  ground- 
reflection  point. 


(3)  The  denendence  of  the  direction  of  the  wave  scattered  into  the 
nrtenna  (relative  to  the  direction  of  the  incident  -..’ave)  on  the 
distance  between  the  cable  and  scatterer  and  on  whether  the  scatterer 


is  inside  or  outside  the  cable. 


~ne  three  effects  (1),  (2)  and  (3)  above  are  all  present  in  the  scattered 
wave,  but  the  n.nthenatical  relationships  are  too  complicated  to  deternine 
exactly  the  relative  role  of  each  of  ther.  Kowever,  it  is  possible  to 
intuitively  Infer  sore  of  this  Infornation  at  least  qualitatively. 

Consider  first  the  "ffcct  ()).  Fros:  the  discussions  in  Fection  5.1, 

•:n  shoul  !  expect  that,  due  to  Tffect  (1),  since  the  scattered  field  arplitudc 


is  rourily  proportional  to  the  illuninating  flel('  anpliturie  at  the  scatterer 
center,  the  scattered  field  amplitude  vill  reach  its  peal:  near  p  -  24  neters. 

This  is  shown  to  be  the  case  for  all  the  outputs  shown  in  Table  5.2,  (a) 
through  (i).  The  peak  is  reached  either  at  P  »  23.5  or  p  =  25  in  every  case. 

This  means  that  (since  computations  are  made  at  1.5  meter  intervals)  the 
true  peak  is  somewhere  near  p  =  24,  i.e.,  when  the  scatterer  center  is  directly 
over  the  cable.  If  F.ffect  (1)  were  the  predominant  agent  In  determining 
the  variation  witn  P,  we  would  also  expect  that  the  1  -  2  dH  assynretry  due 
to  attenuation  (as  discussed  in  Section  5.1)  would  enter  Into  the  results 
'lowever,  tlie  way  in  which  these  features  would  appear  is  quite  complicated. 

There  is  a  contribution  to  any  scattered  field  component  frn-n  the  :t,  y  and 
z-conponents  of  the  illuminating  fields,  so  the  precise  dopennence  of  the 
scattered  field  on  p  would  be  a  composite  of  those  arising  from  ench  illuninating 
field  component.  Since  the  scatterer  is  assumed  to  be  isotropic,  however, 
the  horizontal  and  vertical  components  of  the  scattered  field  would  be  affected 
nrinarily  by  the  horizontal  and  vertical  components  respectively  of  the 
incident  fields;  hence  the  p-variation  of  the  amplitude  of  each  scattered 
field  component  should  contain  some  of  the  features  of  the  corresponding 
component  of  the  incident  field.  The  extent  to  which  this  is  true,  however, 
depends  sensitively  or.  the  direction  angle  of  the  sCBttere<l  field.  Tor  exrrple, 
the  vertical  component  of  the  illuminating  tS-field  will  certainly  give  rise 
to  some  horizontal  components  of  the  senttered  T-fleld  except  possibly  In 
the  case  of  absolutely  pure  backsertter  or  pure  forward  scatter,  which  could 
only  occur  for  n  vary  rare  geometry  of  cable,  scatterer  and  antenna. 

The  affect  (1)  should  cause  rates  of  decay  of  field  strength  sorac'.here 
between  j  and  as  the  scatterer  moves  av;ay  from  the  cable  in  either  direction. 
The  tabulation  below  provides  a  rou'’h  idea  of  wiiat  rate  of  decay  one  might 


exyiect  r'.ui;  to  this  oechanisn.  The  0  rt'5  level  is  ta!:en  to  be  at  0 


in  over"  case. 


0 

3  (approximately., 
nealectl.ne  z> 

?!d-'  '  -20  l°3l0 

i52jdn 

153] dn 

17. h 

5.5 

-22 

-44 

-65 

19.0 

5.0 

-20 

-40 

-60 

20.  b 

3.5 

-17 

-34 

-51 

22.0 

2.0 

-12 

-24 

-36 

23.5 

0,3 

0 

0 

0 

23.0 

1.0 

-6 

-12 

-IS 

2o.5 

2.3 

-14 

-28 

-42 

25,0 

4.0 

-in 

-36 

-.34 

Another  featMre  of  the  variation  of  the  illuninatins  fields  vrith  distance 
fron  the  cable  is  a  difference  between  ^he  illuninatior  mechanise  inside 
the  cabl®  circle  and  that  outside  the  circle.  Inside  the  circle,  the  cable 
-■'boears  concave,  './hile  outside  the  circle,  it  anpears  convex,  "'lus  at  a 
?iven  distance  fron  the  cable,  a  lar'^er  portion  of  the  cable  illuminates 
the  sentterer  than  would  be  the  case  at  that  sane  distance  outside  the  circle, 
'oreover,  inside  the  circle,  the  phases  from  the  illuninatinn,  slots  should 
he  more  nearly  e<lual.  *^his  should  introduce  an  asynmetry  resulting  in  larger 
lllunimtinn  field  strength  at  a  p;ivcn  distance  iron  the  cable  inside  the 
circle  relative  to  that  observed  at  that  sane  distance  outside  the  circle. 

■"he  effect  (2)  should  be  weak  relative  to  (1)  in  deterninin"  the  field 
variation  with,  o  and  is  easy  to  esti'-ate.  Considerin"  the  scatterinn  dirsctl 
into  the  anteran  when  the  scatter  is  inside  the  cable  circle,  the  factor 
— ^  irtrodures  an  i ncrease  in  fiel-'  strer.'^th  as  the  scatterer  moves  awav 
'ron  t'e  ^nble  and  towar'*  the  antenna,  "his  e'fect  is  opposite  to  t''at  due 
to  '.'feet  '1),  whic'  is  a  decrease  in  f^eld  strenth  as  the  scatteror  notes 
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further  inside  the  circle.  The  magnitude  of  this  effect  is  roughly  20  ^SA’’ 

The  nunbers  roue,hly  follov  the  tabulation  below,  neglectin'?  the  effect  of 
the  vertical  coordinate  ?. 


■  '"Si 

-20  1 

23.5 

-27.4 

22.0 

-26.2 

20.5 

-26.2 

10.0 

-25.6 

17.5 

-24.9 

The  variation  betveen  the  values  of  p  for  which  computations  are  nade  is 
only  about  0.6  riC,  The  entire  variation  hctx'een  17.5  ant  20.5  is  only  about 
2.5  (IT. 


As  the  scatterer  recedes  fron  the  cable  when  it  is  outsi'le  the  circle, 
the  effect  is  to  decrease  the  field  strength.  The  tabulation  below  applies 
to  that  case 


‘5.1 


-20  loSio  ''SA 


25,0 

-28.0 

26.3 

-22.5 

2'’.0 

-2?.S: 

decay  of  the 

field  is  about  ,4  or 

1.5  neter  increase  in  0. 

This  ef'cct  is  rou-aly  the  s.nr.e  for  the  eroiind-ref lected  scattered  rie’i, 
since  the  size  of  ir^^  +  r,.,)  is  not  verv  different  Crow  t -.nt  o  '  r^,.  .'.n 

evidences  ’>y  n^ove  tibi’lctions,  it  does  introduce  n  snell  *.*;y-r;ptry  b5f<eer 
of  3  ar  *  outside  t  ie  cable  circle,  '’be  effect  be  to 

rr  !uco  t  ie  rnte  of  docn*  5tre«^*:^  os  tbo  '■c?'tt''ror  'loves  ru  iv  :‘ro* 
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the  caSlo  inair^e  the  circle  and  Co  increase  tne  rate  of  decay  as  it  moves 
wa;'  fro-'  the  cnhle  outside  the  circle.  In  noth  cases,  this  asynmetry  effect 
sho'jld  he  euite  srall  conpared  to  the  effects  (1)  nnd  (3),  the  latter  to 
he  liscusse  !  below. 

’'ionlly  we  consi  ier  the  effect  (3).  If  v/e  thin'.;  of  the  scatterer  as 
be’ nr  illupionted  primarily  by  fields  from  the  nearest  slots,  then  as  the 
sc.ntterer  moves  .nway  fron  tne  cable  inside  the  circle,  the  scattering  into 
the  snhenn.n  direction  hecones  increasingly  more  like  forward  scatter.  As 
it  recedes  from  Che  cable  outside  the  circle,  the  effect  note  nearly  resenbles 
hackscatter.  In  either  case,  when  the  scatterer  is  very  close  to  the  cable, 
it  is  illunirated  priniarily  fron  below  and  hence  the  scattered  wave  Into 
the  antenna  is  at  a  v/ide  angle  from,  the  incident  wave.  Since  hnc-.scatter 
is  likely  to  he  'ciich  stronger  than  forward  scatter,  we  should  expect  that 
this  ■oc'iapisn  will  produce  a  higher  rate  of  decay  as  the  scatterer  novcs 
away  from  the  cable  inside  the  circle  than  that  which  occurs  as  it  recedes 
outside  the  circle,  "hat  asynretry  in  p  around  the  cable  is  exactly  opposite 
to  the  one  cause<l  by  the  effect  (2). 

Tie  res'ilts  on  rate  of  decay  of  field  strength  as  tbe  scatterer  recedes 
fror  the  cable  are  puite  erratic  and  sho-/  no  clear  trends  that  would  Indicate 
whether  ncchnni.sns  (1),  (2),  or  (3)  are  the  prodor.in.nnt  ores.  'T'.e  rates  of 
decay  vary  iron  .snail  values  conrensirate  iflth  nechanisns  (1)  and  (2)  to 
enormous  values  that  appear  to  'jear  no  obvious  relationship  to  the  nechanisns 
we  have  dtscussad  here.  Sonetimes  the  "ic-ld  decays  fron  Its  peak  value  more 
rapidly  inside  the  circle  t‘'a.n  it  does  outside  the  circle  and  sonetires  tbe 
opposite  trend  bole's,  "bis  is  not  surprising  in  view  o'  the  complexity  of 
the  calculations  nade  here,  '■'jich  eiinntity  calculated  is  a  coherent  sum  of 
two  co  ’p'e;:  pii.-ntit  >e3,  the  ffe'd  component  of  the  wave  .scnltored  directly 
into  t'  e  antenna  and  tbnt  oi  t  ’?  scntiemed  wave  reflected  iro  i  the  ground 
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and  then  propagntins  toward  the  antennn.  These  two  complex  nunntitles  are 
of  comparable  magnitude  (as  shown  by  computations  we  have  carried  out  in 
which  these  two  constituents  are  calculated  separately)  and  their  relative 
phase  is  very  sensitive  to  snail  changes  in  pareneter  values,  "ence  the 
resultant  amplitudes  can  range  from  values  near  rero  to  values  nearly  twice 
the  amplitude  of  the  direct  wave  field. 

In  addition  to  this,  the  scattering  process  itself  for  this  frcnuency 
range  is  highly  sensitive  to  angle  and  should  produce  sone  erratic  or 
noise-like"  effects  with  snail  changes  in  geometry. 

At  this  point,  having  discussed  the  variation  with  o  for  fixed  t  and 
z,  we  will  discuss  the  azimuthal  asyrmetry  in  the  results,  hhiere  were  certain 
azimuthal  trends  observed  in  the  illuminating  fields,  presented  in  Table 
5.1  and.  discussed  in  Section  5.1.  For  those  horizontal  field  components 
the  amplitudes  at  i  •  0°  and  o  •  were  nearly  the  same,  os  were  those 
at  h  a  W"  and  h  a  270*,  the  latter  being  larger  than  the  former  for 
and  smaller  than  the  former  tor  “he  z-conponents  were  i^bre  nearly  azi- 

nuthally  uniform.  Small  decreases  in  nmplitude  as  ')  increases  due  to  cahle 
attenuation  wore  also  present  in  all  components  of  the  illuminating  fields. 

All  of  this  was  amenable  to  explanation  through  the  eg-jstions  for  the 
illuninating  fields. 

Since  the  scattered  fields  are  nrogortlonal  to  the  il!”-iir.atlng  f’eldr., 
we  searched  for  the  seme  azimuthal  tren-'s  in  the  former  as  '-ere  o'  served 
in  the  latter.  Some  of  these  trends  are  present  but  not  necessarily  consistep.tly. 
'however,  e:'ccnt  for  some  wild  riuctiintioas  trohahly  ettributahie  to  coherent 
addition  of  two  roughly-e'pial  size  complex  numbers  t*’e  tendency  is  for  the 
a-dlrected  field  conponertr  at  "  ■  'i'*  ■n'-  •  ■  I'K  *  to  he  rou"hly  comparehle 
to  each  ot.her  and  to  the  y-dlroct'’''  corioonents  at  "■  »  TO’  an  '  »  IVO*. 

The  sane  coapcrlsor.  exists  (again  very  roughly)  between  'l,.  at  '  •  CT*  and 
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ar'!  I'^l  at  ''  a  0°  and  nO".  '*o<'ovcr,  the  .t-co'^ponentK  |"_|  appear  to 

exhlhlt  '■ucl'  •'ore  aai-.utaal  as'/'netry  than  their  counterparts  for  the  illunin- 

atinp  fields,  •••itch  '..ere  nearly  svrrotrica! .  "here  is  also  noro  azinuthnl 

asynnecry  in  all  conponents  of  the  scattered  fields  than  i.’as  orosent  in  their 

counterparts  in  the  illuninatinj  fields. 

A  nossible  partial  reason  for  this  additional  asynnctry  in  all  the  conponents 

is  the  f.nct  that,  -at  d  •  0°  end  0  =  1?0“,  the  scattered  field  conponents 

contain  Irryo  contributions  fron  E  and  "  of  the  incident  fields  and 
X  '  z 

the  scattered  field  conponents  contain  Inrse  contributions  fron  incident 

and  '  ,  Also,  at  c  „  <50“  and  *  =  270“,  scattered  1  has  contributions 
z  y 

fror  incident  T,.  and  contributions,  •'htis  asynretrles  in  ,1,,  and  il..  in 

tic  illurinntino  field  wil'  find  t'leir  ■.•ay  Into  ^  jn  the  scattered  field 

z 

and  will  also  introduce  greater  asynretrles  into  scattered  ^  and  !!  . 

•:  y 

or  're  nest  exanine  the  differences  between  the  results  for  the  radial 
"nr-  (z  »  1,  upri;‘'t  .acotteror,  Tihle  5.2,  (a),  (b),  (cl),  the  radial  cr.Twl 
<r  «  '  .1’ ,  rrono  scntterer,  "able  5.2,  (d),  (e),  ff))  end  the  radial  wall; 
on  stilts  (z  ■  2,  uprieht  scatterer.  Table  5.2,  (p),  (h),  (1)). 

Aside  fron  the  trends  already  cormented  upon,  which  tend  to  be  connor. 
to  all  throe  cases,  there  arc  sone  general  conparative  tendencies  observed 
which  nnnesr  to  be  consistent  with  intuitive  expectations. 

"iryt  conpari.i'’  tt'e  radial  wrlh  on  stilts  (z  »  21  with  the  radial  wnl': 

’  »  1),  the  correspond iny  field  conponent  annlltude  for  the  forner  tends 
to  hr  snnl'or  chan  that  for  the  latter.  This  is  sonctines  strikingly  true 
an  '  •'onoti-’'’S  the  reverse  holds,  so  one  cannot  say  that  trend  is  absolutely 
cons'rt''nl.  ’’ut  it  is  true  nore  often  tlian  not  and  tends  to  be  accentuate' 

';  •  pro:ti  itv  to  f-r  ca'ile,  i,c.,  it  is  nore  likely  to  be  true  for  the  pe.n!: 
values  I'-.op  for  the  vanes  '•  -  5  -eters  fro-,  the  c'ihle.  An  artplanation  is 
test  t  le  effect  o^  the  hoie'*t  (the  v.olne  oi  z(  is  ruch  nore  oronojincod  verx- 
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near  the  cable,  bein;  comparable  to  the  horiaontal  distance  from  the  cable, 
and  becomes  nearly  neslly.ible  6  meters  fro'r  the  cable,  vhere  the  heie''t  in 
both  cases  is  snail  compared  with  the  horizontal  distance  from  the  cable. 

T!ie  sane  observations  can  be  made  concerning  comparison  bct'.een  the 
radial  crao!  (z.  ■  0.25)  and  the  radial  i.-elh  (z.  ■  1).  Tor  most  of  the  points 
at  vhich  connotations  './ere  made,  the  radial  cravl  <’ives  larger  values  tlnn 
the  radialv  alk.  This  is  particularly  true  of  the  pea!:  values,  which  are 
nearly  always  significantly  larger  for  the  crawling  case  than  for  the  val!:ing 
case.  The  field  magnitudes  3-6  meters  away  from  the  cable,  however,  are 
not  always  larger  for  the  crawling  case  and  are  sometimes  sorowhat  smaller. 

To  some  extent,  this  can  be  attributed  to  the  smaller  sensitivity  of  the 
amplitude  in  the  crawling  case  to  the  height  z.  !:lnce  z  for  the  crawler 
is  one-‘’ourth  of  its  value  for  the  walher,  it  beconer,  ne"ligi')le  compared 
to  the  horizontal  distance  from  the  cable  when  the  latter  is  only  nhort  2 
"■•oters.  hence  over  much  of  the  range  of  o,  the  comt-arptive  nmplitu-'es  in 
the  tvo  cases  have  very  little  to  do  with  the  concaratlve  values  of  the  he’cht 
ithin  two  meters  of  the  cable,  however,  there  j_s  n  signi'^iennt  dependence 
o'  the  amplitude  on  z,  and  than  helps  to  ezplrie  the  fnct  that  per:  va’uoE 
of  c"pHtude  are  nearly  always  significant';/  higher  for  the  crawl  in"  case. 

'side  from  the  comparison  involving  the  height,  nlludei  to  aitove,  the 
remnindor  of  the  comparison  het>/ccn  cre’/ling  and  wa'hir.g  cases  relal'-s  to 
tlie  difference  in  the  angular  orientat  ion  of  f'C  sentterer  in  the  two  cases. 

In  the  vT’!:ing  case,  the  scatter  in"  tends  to  be  "bronrside"  '.'hereas  In  the 
crawlin"  case  it  tends  to  ^e  in  the  snhoroiJ's  lon"itt!dinal  -'irectior.  "he 
scniteria;;  in  the  latter  case  tends  to  be  weaher  than  that  in  the  former 
case;  ''.erce  for  so"?  sentterer  positions,  tp,.  feet  t''at  the  ill  nr  iraeir" 
fie' i  is  stronger  in  t'’e  cr.a'.lir"  case  f./oceute  th"  corner  of  the  rra -.'er 
is  'o  /er  and  tlieroforc  closer  to  the  c- t 


'■"n  t  vt  of  the  ''.-.r  ’r  nt  t' r 


sn-;''  horizontal  Usr?nco  fron  the  cable)  is  offset  by  the  fnct  that  the  scatterin; 
is  «;tron?or  for  the  walker  due  to  its  orientation.  That  in  more  li!:el>-  to 
be  true  far  frnn  the  cable,  where  the  heieht  difference  has  little  effect, 
uhil?  th'-‘  '’ei7bc  is  likely  Co  play  the  dorieanC  role  nenr  t^'e  cable, 

•.b'-'  ori?nt‘>tioa  should  be  ror‘'  i'^'ortenf  further  ^ro*^  the  cable,  ’.'hich  ai^'-t 
e::plrin  •■h’  the  pcact'‘red  field  is  soretines  significantly  larger  for  the 
\ir]  r>T  ospecinMy  at  positions  far  fror  the  cable. 

*s  a  tirrl  observation,  note  that  the  horizontal  and  vertical  components 
of  the  ncettered  field  arpHtudes  are,  on  tae  average,  not  radic.nlly  different, 
“•Ithoua'i  ^h'';.'  do  exhibit  scat'  erratic  behavior  (e.c*. ,  in  some  cases  there 
are  !ar  ,c  differences  in  the  noak  values  for  different  components,  v^here 
vt-rticil  is  soT^tiros  larger  and  sometimes  smaller  than  the  horiaontel  com- 
nofirnts).  'h-*  •'robable  causes  of  t^'e  erratic  or  noise-like  nature  of  some 
o^  » be  re^iultf  have  already  been  liscusscd,  allurin';  to  the  sensitivity  o'^ 
the  scatterin'?  process  to  small  chan'^os  in  geometry  and  coherent  addition 
ronai'':  Timbers  vith  sensitive  relative  pltasc  as  the  major  rechanisns 


CTjsirm  t'lls  behavior. 


5.2.3  Total  Fields 


The  computed  results  for  the  comoonents  of  the  total  field  at  the  antenna 
(fields  without  Che  scetterer  plus  those  due  to  the  presence  of  the  scatterer) 
are  shown  in  Table  5.3. 

Discussion  of  the  fields  due  to  the  presence  of  the  scatterer  wns  covered 
in  Section  5.2.2.  The  computed  field  conoonents  at  the  antenna  in  the  absence 
of  the  scatterer  are  ?iven  by  (5.1)  and  are  -112  IB,  -110  dB  and  -Idf  dB 
for  the  X,  y  and  z-electric  field  amplitudes  respectively. 

These  values  were  determined  by  a  very  accurate  numerical  technique. 

That  technique  and  some  approximate  calculations  of  the  field  components 
are  discussed  in  Appendix  I,  where  the  fact  that  the  horizontal  components 
at  the  antenna  is  much  larger  than  the  vertical  component  is  explained. 

"he  remaining  discussion  in  this  section  will  be  confined  to  comments  on 
the  parameter  regimes  in  which  the  field  component  amplitudes  given  by  (5.1) 
for  the  case  where  the  scatterer  is  not  present  nre  comparable  in  mo’.nltude 
to  those  given  in  Table  5.2  for  the  contributions  due  to  the  presence  of 
the  scatterer.  In  those  regimes  there  is  a  significant  effect  attained  through 
coherent  addition  of  the  antenna  fields  in  the  absence  of  the  scatterer  and 
those  due  to  the  scatterer 's  presence.  In  many  regimes  the  scatterer-froe 
conoonents  given  by  (5.1)  overwhelm  the  scattered  fields  .and  the  results 
show  the  scatterer  to  be  undetect.able.  Tor  other  regimes,  fie  components 
given  hv  (5.1)  are  overv;helmed  by  the  scattered  fields  and  the  values  in 
T.able  5.3  are  almost  ex.ictly  the  same  as  the  corresponding  values  in  Table 
5.2. 

To  initiate  the  discussion,  we  sommarize  the  results  in  the  tahuloiior. 
below  (obtained  from  Tables  3,1  and  5.2  with  the  aid  of  Tq.  f'l.D)  which 
indicates  the  points  of  detectabi' ity  nf  t‘’e  scatterer,  denoted  :iy  and 


characteri^c-i  the  fact  that  the  presence  of  the  scatterer  changes  the 
field  cornonent  at  the  antenna  by  at  least  4dp;  an  incret’ent  below  4  dB  should 
;>e  considered  as  ’’noise.**  As  a  subclass  of  these  points,  we  also  indicate 
points  where  the  scatterer  contribution  to  the  field  overwhelms  the  field 
cornonent  that  would  be  seen  in  the  absence  of  the  scatterer,  i.e.,  introduces 
an  incrofTont  of  ar  least  *  10  db;  these  points  are  denoted  by  '^A.  Points 
where  there  nr<'  no  perceptible  changes  in  the  antenna  sipnal  due  to  the 
presence  of  the  scatrerer,  i.e.,  -here  the  scatterer  Is  "not  detectable” 

''i.r*..  ■_  ](’'  iP  change)  are  denoted  by  'ih.  '“inally,  those  points  where  there 
IS  ^  smell  nut  harel  v  nercentat'le  change  in  the  sipnal  level  due  to  the  scatter 
“r's  pres#»nce  "barely  detectahle"  scatterer)  are  denoted  bv  BH;  between 

1  '.?•  aid  4  dl  change  due  to  presence  of  scatterer). 
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Thf!  issue  of  detectnbllity  of  the  intruder  is  that  o:  the  conparison 
b?tween  the  anplitude  of  field  component  at  tne  .intenna  v.’ben  the  senttorer 
is  present  relntlve  to  the  anplitude  of  t*'at  sane  field  conponent  vhen 
the  sertterer  is  not  present.  In  some  rises  the  total  field  a"*plit»de 
i.s  reduced  rather  than  increased  due  to  the  presence  of  the  •: erer .  The 
re'^uction  will  occur  if  the  amplitudes  vitii  and  without  the  sentterer  are 


of  conp^rable  T^a^’nitude  and  the  relative  phase  an»le  of  these  two  quantities 
is  between  10®  and  270®,  while  an  enhanceirent  occurs  when  that  angle  is 
between  -00^  and  -^90®. 

Tr  the  existing  systen,  the  antenna  is  vertical  and  hence  the  z-component 
of  the  cJectric  field  is  the  onlv  ronnonent  of  interest.  However,  there 
IS  no  r''3Son  fin  principled  that  horizontal  conponents  could  not  be  used 
for  detection  with  an  antenna  that  respond  to  those  coirponents.  Therefore, 

It  IS  worthwhile  as  a  part  of  this  study  to  consider  the  horizontal  conpon- 
epis  of  the  field  in  addressinp  the  detectability  issue. 

”he  overall  conclusion  that  wouH  follow  fron  the  results  shown  here  is 
that  toe  dtcctahilitv  of  the  intruder  with  horizontal  fields  in  all  three 
of  the  cases  considered  t«  nonexistent  or  very  ooor.  In  a  few  cases  involving 
the  radial  valk  or  radial  crawl  and  one  case  ’./ith  radial  ua^k  on  stiles, 
ch?  scatterer  is  brrely  detectable  at  a  point  n^ar  the  cable,  but  in  nearly 
<111  ca^es  it  is  not  detectable. 

If  the  vertical  field  coinponent  is  used  for  detection,  as  in  the  actual 
systen,  the  conclusions  that  would  follow  fron  the  results  is  that  the 
intru<ler  ls  easilv  detectable  for  the  radial  walk  and  radial  cases,  but 
sonew‘'at  less  detectable  or  undetectable  for  the  ”r<adial  walk  on  stilts” 

"ase . 

Since  /e  have  experirental  results  for  conparison  only  for  the  vertical 
^'lectric  field  conponent ,  no  further  connents  on  this  topic  v/ill  be  made 
concernin'’  the  Horizontal  field  conponents.  The  oxoerirental  results  for 
rpp  vertical  antenna  tisial'y  sho'.*  values  of  the  mnlitude  o?"  the  scattered 
si'^nal  at  the  antenna  sonewnerc  b'^tveen  9  and  10  f*o  above  the  signal  in 
th''  absence  t*  t^'o  scetteror.  Oyr  analytical  results  usually  show  a  larger 
inrrenc'^til  s’^nal  clue  to  the  scatterer  than  in^-'calo'  the  exoori'nental 
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results,  i.e.,  for  the  radial  ua'h,  the  peak  scattered  sienal  values  are 
between  17  and  31  d?  higher  than  the  antenna  signal  in  the  absence  of  the 
scatterer.  For  the  radial  crawl,  they  are  between  13  and  48  dP.  higher 
and  for  the  radial  walk  on  stilts,  they  are  between  10  and  22  c'P  higher. 

'.■hen  a  scatterer  is  a  few  neters  away  fron  where  the  peak  values  occur, 
the  increnental  signal  due  to  the  scatterer  is  nuch  more  nearly  co.mparablo 
to  that  without  the  scatterer. 

However,  although  the  qualitative  hehavior  of  the  signals  due  to  the 
scatterer  .appears  to  be  roughly  as  expected,  the  most  prorinent  quantita¬ 
tive  descrepancy  appears  to  be  excessively  high  values  of  the  scattered 
signal  relative  to  the  signal  In  the  absence  of  the  scatterer.  This  is 
either  because 

'1  The  calculated  field  amplitude  at  the  antenna  in  the  absence 
of  the  scatterer  is  too  low,  or 

(2)  'be  calculated  field  amplitude  dje  to  the  scatterer  is  too  nign, 

txplanatien  (1)  seems  less  likely  than  explanation  (2).  '■'he  'ield 

amplitude  at  the  antenna  in  the  scat'erer's  ab.sencc  was  a  stralghtforierd 
and  simple  auantity  to  c.ilculate  and  presented  no  major  conputntional  pro¬ 
blems.  fic  calculation  Is  discussed  in  Appendix  I,  where  opnroximate  methods 
are  shown  nnd  where  hlgaly  accurate  numerical  integration  methods  that 
were  actually  used  to  obtain  the  results  are  alluded  to.  In  doing  this 
computation,  the  effect  of  the  ground  reflections  was  found  not  to  be  as 
significant  as  in  the  case  of  the  fields  illiirlnattng  the  scatterer,  i.e., 
the  fields  at  the  antenna  with  or  without  ground  reflections  are  cor.parabla 
in  magnitude.  The  effect  ground  reflections  could  have  liad  woul'*  be  in 
the  direction  of  reduction  of  the  antenn.a  field  due  to  destructive  Interference 
between  'inset  and  ground-reflected  si-nuls,  or  in  the  direction  o'  an 
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'?nhnncr‘'^'nr.  'n  >  fp*  due  to  constructive  Interference.  Therefore, 

althnut^n  wcokne^^s  in  iccountinp,  for  r.round  reflections  is  a  possible  cause 
of  inocc  Lr.jc in  this  ro'^ult,  the  ?ffer^  of  inprovine  that  part  of  the 
rnlculati^n  •ni-'.hr  Keen  in  the  -.^ronr  direction,  i.e.,  the  direction 

oT  r':''iucti'^n  rnthor  than  enh.inc'^'^ent ,  or  at  best  an  onhancenont  of  less 
’’hnn  'I  !’■. 

*' ';~il  arat  ion  is  more  likely.  It  appears  that  the  spheroid  used 

-iodr>i  the  intruder,  t-ith  the  oaraneter  values  */e  nave  c^'csen,  is  too 
’.ood  a  scntt‘=*rer.  Once  r  decision  had  been  mado  on  the  choice  of  parar?eter 
values  and  rhe  for'iot  for  the  comutations,  there  was  inadequate  tine  or 
resources  to  na  -e  chan'cs  in  these  vah’es  and  thereby  determine  enpirically 
which  values  vould  provide  “.he  best  fit  to  experinental  results.  There- 
fore  the  cor-)«-.tations  \  pre  all  made  with  these  same  values.  The  values 
c  h os 0 r ! 

=  '.npth  of  spheroid  *  2  meters 

"  a  ''adius  of  epaeroid  *  0.25  meters 

*  Permittivifv  oi  spheroid  •  ,4080  (10  )  farads/rcter  *  46,2 

c  s  Conductivit?  of  spheroid  *  ,502  »hos/meter 

''' c^  opo  of  those  parameters  could  have  been  chosen  to  have  a  gnaller 
‘nc''  •••m-jld  hav#»  produced  a  smaller  scattered  fi^'ld.  The  values 

ictuallv  lar'-er  than  those  corresponding  to  a  nan  of  overspe 

«;  o 

S’^c  is  -^pu'valent  to  6,5m  feet,  extremely  tall  for  a  man  and  3'’^  is 

enuivaiert  r<  ’.64  f?et,  extremely  wif'e  for  a  man). 

’  chosen  valics  o’  c  and  0  were  bosc'^  on  consultations  \/ith  ’’rofessor 
s  s 

'nrher  and  are  hared  on  composite  values  from  constituents  of  the  human 
body,  “her?  ’  rrob.n'nU  vide  variability  in  the  optimal  choices  of  tne 
'■onst  i  t ’It  ive  ^nnneter's  for  this  scatterin;’  model,  Tedurtion  in  either 
or  hot’  -'o-il.  hivrt  ’■oduc'^d  the  .scattered  fields,  '’^e  other  un''nov.-n 


is  the  corresponriance  between  the  actual  height  and  girth  of  a  hu'ian  fra"te 
and  the  proper  height  and  radius  of  a  spheroidal  scatterer  used  to  nodel 


a  huinan  frame. 

The  summarize  the  point  of  the  discussion  above,  some  degree  of  "cut 

nd  try"  choice  of  the  scatterer  parameters  L  ,  R  ,  s  and  C  night  be 

s  s  s  s 

necessary  to  natch  the  calculated  received  signal  at  the  antenna  to  experi¬ 
mental  results  before  the  codel  could  be  used  to  predict  results  o‘  changes 
ip  system  design  parameters.  Other  parameter  values  used  in  the  computations 
are  much  easier  to  choose,  because  they  are  usually  fixed  by  the  geometry 
of  the  system. 

"here  are  many  approximations  uaad in  this  model,  and  there  are  noise- 
lihe  fli:ctuations  in  somu  of  the  computations  due  to  the  sensitivity  of 
the  scattering  model  to  the  geometry  and  to  coherent  additior  of  complex 
numbers.  The  relative  phases  of  the  contributions  might  he  subject  to 
numerical  errors  and  the  results  arc  extremely  sensitive  to  these  relative 
phases,  "here  are  also  fluctuations  In  the  experimental  results  (e.g., 
fluctuations  as  the  ancle  f  chor.ges  in  the  "circumferential  wal''."  results), 
so  it  is  no  surprise  that  such  fluctu.ntions  occur  in  the  analytical  results, 
fiouever,  general  trends,  such  as  the  peol-.ing  of  the  signal  '/'-.en  the  scatterer 
is  near  the  cable  and.  decaying  as  it  recedes  from  the  cable,  ir?  present 
in  the  results,  it  is  concluded  that  this  analytical  model  and  nrsociated 
computer  program  could  be  used  to  predict  the  effects  on  performance  of 
changes  in  environmental  parameters  (i.e.,  different  types  of  soil  and  design 
parameters  e.g, ,  different  cable  dimens’ons  or  antenna  position.s  or  orientations) 
or  intruder  parameters  (i.e.,  height,  girth  and  constitutive  paroneters 
of  scatterer).  In  its  present  form,  the  computer  program  contain'  wide  flexi- 
h'lity  in  parameter  value  choices  snd  La  rot  excessively  computer-tire 
i ntens’ ve. 
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'uperpoi^icion  of  the  field  contributions  fron  slots  nrounc'  the  cable 


is  accomplished  through  the  evaluation  of  the  integrals  given  by  (2,17-a 
h,  c''  with  cho  aid  of  (2.1^)*  In  this  appendix  we  present  some  of  the 
detail*^  concerning  approximations  tor  this  integral  (and  for  the  field 
components)  when  rno  observation  point  is  close  to  the  antenna  (Section 
T-A)  and  when  it  is  clos'^  to  the  cable  (Section  T-^).  In  Section  I-S, 

Che  nurorical  soUtcion  for  aribtrnry  location  of  the  observation  point 
ir  briefly  discjssed*  In  Section  ve  discuss  some  methods  such  es 

stationary  phase  vhir>  can  be  nsed  to  ovaluft?  the  integrals  appro::imately 
for  nrbitraf’'  observation  point  locations* 

1-A  Ipproxinaticns  for  *)bservn*  ion  Joints  bear  the  Antenna 
'ron  the  -e‘'inition  of  ®  Selo*  Tqs  (2,l-o,  b) 
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■^■he  I'-nroxi-.jt  10'  that  aoniics  to  observation  points  near  the  antenna 
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wiich  is  true  for  oil  case?  of  interest  in  this  study. 

^  n 

Ne^lectino  the  pecond  term  in  (1-5  ir  the  omolitude  (»“  )  and 

retnining  it  in  the  phase,  we  cnn  write  the  inte'jrels  fron  (2.17-a, 
5,  c)  with  the  nii  of  the  approxiration  (7.13)  cs  ^ollo's: 
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In  the  case  where  the  ohservntiot!  potnt  ie  at  the  e:MCt  rente'  o‘  the 
cable  confl'’.i'ration,  i.e.,  o  »  C,  the  Ilesacl  'unctionr.  ether  than  Jq'’*'’' 
vanish  an'i  the  integral  (1-5)  is 
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""his  is  t  ie  case  that  annlies  wJier  one  calculates  the  ficH  at  the  antorm 
when  the  antenna  is  nlacei'  at  the  exact  center  of  the  crb.'e  confi"iiralion. 

The  ersicst  wav  to  calcjlnto  the  fiaid  components  for  this  case  i? 
to  he'in  ’■>  usin';  (?.’rh-a,  b,  c'  an-I  (2.11-r,  b,  c)  to  calculate  the  conpon- 
erts  o;  the  fields  due  te  a  sin'le  slet  at  anrle  n'.  'e  re'rite  these 
eei'ations  in  ter'S  of  cot  in  ter's  ol  an'  r"',  thus  r'?-ie,vir  t''e 
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•Tonaront  f^cpenJence  on  cho  observation  point**?  azimuth  an'^le  rt,  Tn  this 
particular  case  the  fields  are  indepeocent  of  0  ftruc  for  no  other  case) 
nnc  the  formulation  in  terns  of  ,  useful  in  all  other  cases  should  be 
r^viiced  hy  a  fornulation  in  which  all  quantities  are  expressed  directly 
>r  rer*^*^  of  '*j'  ratlier  than  in  ternr  of  *'  and  v*. 

T  ^  Hosir<*d  forns  of  (2.10-a,  h,  cl  and  '2.11-a,  b,  c),  usipo  the 
■*n"*ro'<  i"’at  f  or 
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''innlic'!  'jv  M-2),  '1-1)  and  the  ’.no»n  value  of  A)  can  be  approximated 
bv  the  oxnrcsslons 
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Integrpting  b)  fror  *’  •  0  to  •  2t  and  involting  the  aoprox- 

imatioi  (2.1?)  in  evaluating  these  integrals,  vie  obtain  (approximately) 
the  fom  (1-7)  for  all  the  integrals,  leading  to  the  result 
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7tg  results  (T-lO-a,  b)  provide  ronvonient  ipproxinat ' nn'^  th'^r  enn 
■)e  jser  for  connnrison  with  the  nuTiericil  ly  ev3luat3<*  fieJI  corponont? 
ob  tained  for  the  cise  where  the  sentterer  is  absent  and  ", round  reflections 
are  not  accounted  for, 

‘^’be.«e  results  shov  th«  predonin-mt  electric  riolri  conponert”  in  Lie 
horizontal  lirectior  and  oredc^^inrnt  »“nfjnetic  Me\d  CDT-'nonent  in  tl^" 
T-ijrectior.  The  l-'tter  re?  ill  eonmr^  to  he  Hue  to  t'r  fret  th-t  the  •st* 
netlc  fizl"'  in  th*  cchle  slots  is  a^l.-uthell y  <’irccto:’  '’.e,, 
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results  obtaised  fro't  the  rigorous  conputations  usin”  a  highly 
acente  n  j-erical  ion  teclmipue  (See  Section  T-S)  arc  p,iven  by 
5,1  ?n  i  are  (  roiin'led  off  to  the  nearest  inte-er) 
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^iico  r.r«  inteprrtior  carried  out  in  this  appendix  for  the  case  0  »  0 
\r,  evact,  ir  ir?  'mr<ily  surprisin';  that  the  a^reenent  between  (1-12-8,  b) 

'•’0'!  (l-l-T,  b)  is  so  'locd^.  fowever,  it  is  still  useful  to  establish  that 
goi-Artncnt  bccnusc  it  serves  to  stren^hten  cQnfi,tenco  in  the  nunerical  lethod 
lo  e''nluat':  the  into'^rals  in  the  '•.eneral  cnr»c. 

I-  ionrox  i'^ation  Cor  Observation  Joints  *^car  the  Cable 

The  ir^’pnrtsnt  iren  of  int^'rest  for  the  scatterer  location  in  this 
sC'id]'  is  re'^ior  near  the  rrble.  To  ^’etenine  the  fields  »t  the  scatterer 
’orntlcn,  t^?rrfore,  it  is  useful  to  ’ocus  on  the  approximations  valid 
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which  is  the  radial  distance  fron  the  cable,  positive  inside  the  cable 
and  nejetive  outside  the  cable. 


’'e-writin'’  (l-l)  in  terrs  of  A,  ve  obtain 
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'JnfcT-timately  it  is  difficult  to  use  thi?  aporoxi"ation  effectively 
unless  Z'  is  Inree  enough  so  that  the  second  ten  inside  the  rnnlcnl  in 
fl-l-il  is  snail  conpare'*  with  tie  first  tern,  is  the  ancle  t'  increases, 
this  condition  becones  increisin",! y  vnlid. 

fhere  are  three  regions  of  interest  '.’ith  resnoct  to  tl>e  depenler.ce 
of  ■'  on  the  anjle  'I'.  To  discuss  this,  ue  define  the  second  ten  jesi'e 
the  ndica’  in  (1-12)  as 
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"ith  inranot^r  value  3S5i ’nnent  •  ir  our  connotations 
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ri:  sin^  ^  •<■  (z  -  b')^  (I-17-a) 

^2:  sln^  ^  ~  >L? i  is  as  "ivea  by  (I-IA)  (I-17-b) 

■■3;  sin^  ^  »  rv.‘;  R  is  2|l|/l  jsin  (I-l7-c) 

I - 2  -r 

’otia'’,  tbat,  ir  3erRnl,  if  x  <<  1,  then  1  +  x  -  ^  +  y  >  it  If 

ovUlent  that  ^  valae  of  x  equal  to  y  will  produce  an  error  of  ,  or  roushl 
3  nercont,  If  \/l  +  x‘  is  approximated  by  I,  Tlie  rocions  arc  chosen  accorrl- 
ln®ly,  i.c.: 

Rl:  jsin  <  y 


n  I  o'  i 

02:  ^  <"  |Str.  -^j  <  57L 

■’3:  jsin  >  5?t, 

ohc  tch'jlatioa  helo-.’  Indicates,  fcr  dif 'erent  values  of  the  maximu"! 
•in^le  for  the  resion  21  and  the  mirimun  an",'e  O'  for  tlie  region  ?3.  Thes 
ir’ler  are  ineicnted  hy  Oj ,  and  '.I 


respectively 


8 


- 

radians 

de'>raes 

(f'^)  radians  (Z 

p  degrees 

-6 

.0452 

2.590 

1.20 

5.'.r 

-4 

.0315 

2.127 

.012 

45.5 

-- 

.mo 

1.031 

.454 

26.0 

-1 

.0012 

0.642 

.230 

15.0 

0 

.00.30 

0.453 

.200 

11.5 

1 

.0120 

0.533 

.300 

17.2 

2 

.0103 

1.123 

.490 

2'?.0 

4 

.0373 

2.154 

.820 

47.0 

6 

.0534 

3.346 

1.46 

33.7 

Tho  purpose  of  the 

above  discussio*' 

is  to  define  regions 

along  the 

cab  I 

e  here  the  integration  over  that  portion  of  the  cable  can 

be  easily 

anrroxlmatert .  The  approxiiate  int'j'rnl 

over  81,  froi  (2.17-a) 

.  ^2.13) 

an^' 

n-17-a).  i.s 

•  c 

-  •  I" 

-  > 

eJ'^n  - 

-(a  +  jkf'  )£' 
ca 
e 

;  (o  e  - 

+  e 

2iT)  )#2  r 

-  0 

(i-i  ) 

•■’ifcre 

1^7  7 

■?L  -  I  -  (*  -  b’)- 

V}'}  iTnror'rwte  irtesrrl  ov?r  T?,  tron  ^2, IS)  anH  'I-i7-c), 

is 


r-o 


"nfortunate]  y,  the  inCe’nls  in  (I-l'')  cannot  be  evelurted  exactly,  nor 
rrn  the  inte'-rnl  that  ■..•oul"'  arise  for  region  P2,  where  no  nnpro:  ir.ntions 
li':e  tho'-e  ese''  in  (I-l'-'  or  O-l'’)  can  be  nsc-i.  icnco,  nnlers  wo  can 
validly  nsr.iiae  that  the  inteyral  yiven  in  il-lC)  is  the  over'.helrina  contri¬ 
butor  to  the  irce'’ral  over  the  entire  cable,  the  aoeroxi-i? t i ons  due  to  the 
close  arexioity  of  t he  observation  point  to  the  coble  do  not  help  to  evaluate 
the  integral.  The  assurption  that  the  inte"rel  in  (I-IT)  is  a  pood  approxi- 
aation  to  the  inteernl  over  the  entire  cable  is  piven  credence  by  sone  of  the 
aaproxia.atc  results  in  ^action  5  (i.e.,  those  based  or  Table  5.1  and  F.qs. 
(5.3-3,  o))  wtiich  seen  to  indicate  that  the  overv/helnin-  contribution 

to  the  11  lur-laat ing  fields  cones  fron  the  slots  nearest  the  observation 

point  and  thnt  for  each  component  the  field  annlitiides  seen  to  vary  rouplily 

1  2  2 

as  a  pnrticul.ar  power  of  4,  where  o  »  +  (x  -  b')  . 

I-C  The  fleneral  Case 

"he  actual  evaluation  of  the  inteprals  inplencnted  for  nachine 
co“putation  ns  acconal Ishe  '  thronph  a  highly  eccurate  ninerical  integra¬ 
tion  tec  'hieuo.  d  t'  ulftii  order  Gaussian  rule  (^onte  and  de"’oor,  ""leaentary 
niinericnl  analysis  -  an  nlopritbnic  anprocch,"  "cCraw-'Ii  1 1 ,  1 1E0;  Pages 
311-31.i,  325-32b)  v/as  used,  "^e  first  12  coefficients  in  a  series  are 
deterrined  in  order  to  evaluate  the  teros  of  a  set  of  12^  order  orthoponcl 
nolynorialn.  'sinv  12f''.'  to  2f0f  points,  the  nethod  provides  accuracy  within 

3.1  po'cent.  i’he'^e  coapitat  ions  i  ere  carried  out  •..•ithoul  significant  expend.i- 

fnre  o'  rurnirp  tine  and  tables  of  the  vclues  of  I  for  reouirer  values 

n 
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of  0  and  I  were  generated.  This  task  was  sitpliflod  by  the  fact  fiat  the 
integrands  of  depend  only  on  G'  -  b*  -  b  and  are  Independent  of  n,  the 
azlruthal  angle  of  the  observation  point. 

T-h  General  Amroxinations  for  the  Interrnl  1^ 

Tie  integrals  given  by  (3.17-a),  (2.17-b)  or  (2.17-c)  (appro?:i  ’..ite’y 
eguivelent  by  virtue  of  (2.13))  are  all  of  the  generic  form 


-  p  dz  ;(x) 


(I-2n-.a) 


where 


(•dx 

'  V 


f(k) 


(I-20-b) 


:<  -  O' 

l(x)  •  -  2  p  cos  X  *  (,t 

1 


f(xl 


V)' 


(^(x)r 


;  n  ■  1,  2  or  3 


V(x)  >  kCCx)  -  k^^di.*) 


e”'^"'  f(x) 


-J.r) 

*(x)  .  jf(x)  -  ri 


'^'.c  forn  (I-2fi-a)  i.’ill  be  i.'set’  in  the  disrussion  of  the  stationary 
n'-ase  retiiot'.  The  ‘'orr  (I-2'l-b)  n.ll  be  used  in  discussion  of  a  technique 
that  is  an  alternative  to  st.'’tionarv  g’  ase  In  cases  vbere  r.o  stationary 
phase  ''oiats  exist  or  the  latter  aethod  ta  not  easily  applicable  for  other 


reasons. 


✓  /  ^ 


w(*  r'now  that  ^  1,  which  is  a  contrar^iction  of  (1-30)  or  (1-30)*; 
therefore  the  positive  sign  in  (J-26)  is  outlawed  by  the  requirement  that 
cos  £  1 .  which  we  have  shown  to  ho  inconnetible  with  that  case. 

consider  the  noparive  sip.n  in  (T-26),  which  allows  cos  to  be 
neoativo  ir  positive,  '"he  required  condition  is 


[cos  XqI  <^1 


(T-33) 


or  equivalently 


cos  £  1 


(1-33)* 


'  aich,  hasod  on  n-26),  has  the  forn 
d')  ^  2 


/  d2  .,2  /  2  \  f  2 

^  .  1  (pz  .  ,,  •_  ’£^ .  1  -  1  <  1 

'  0  2  I  ca  0  1  2  \  ca  — 

If  )  V  0 


'ron  (1-27)',  we  can  write 


P  -  T~ 

ca 


,  where  3  1 

''sin?,  fl-3j)  in  (1-33)",  we  have  the  condition 


(1-33)” 


(1-34) 


:2,'2 

ca 


‘  -  l)(t^  -  1)  -  23.'  A-  -  1  - 

ca  ca>/  '  ca 


) 


•5  a?  A  M2  I  2 

:3'(2!<‘-  -  1)  -  2!;  .k*^  -  1  S  .'C  -  1  -  <  0 

cn  ca  \  c£  ,  ca  - 


'>2 


(1-35) 


sin?  the  nsslijned  value  o'  »  A. 7,  we  con  re-vxite  (T-35)  in  the 


-  2.13S  I  3-  -  1  -  1.7  £  0 


(T-36) 
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A  study  of  the  Inequality  (1-35)'  shows  that  there  Is  a  liaiced  ranqe  of 
values  of  S  (and  hence  of  0)  where  stationary  phase  nethods  can  he  used. 

Since  stationary  nhase  points  do  not  exist  for  all  values  of  P,  it  vas  decided 
not  to  use  that  method  to  evaluate  the  integrals  An  alternative  approxi- 
nation  method  for  integrals  of  the  type  (I-20-a)  or  (I-20-h)  will  be  discussed 
in  what  follows. 

■•,'o'i'  we  consider  an  asvnntotic  method,  which  was  used  in  the  original 
project  and  is  given  to  first  order  in  Reference  1,  Section  4,  Page  4-3S, 

"Tiis  method  applies  to  cases  not  easily  amenable  to  stationary  phase  methods 
and  is  based  on  partial  integration.  To  derive  the  method,  we  use  the  form 
fT-20-bl  and  integrate  successively  by  parts,  i.e.. 


’‘i  *1 


After  two  more  steps  of  Integration  Sy  parts,  we  obtain 


-inisii  , Il-J— ;  .i: 

[■.’(xu'  -■  ' 

•  -Jl] 


4"(x)  _ 

Ih'Cx)]^ 


This  method  can  be  applied  to  the  integrals  ns  long  ar.  'ocs  not 


vanish  at  the  end-aoints.  ’lote  that 


and  hence  fron  (1-21) 


'tVx) 


•  sin  » 


-  k  )  -  a) 
ca  ’ 


(1-39) 


Tf  the  ond-points  are:  Xj^  “0,  ^2  *  then 

-  'S’CXj^)  .  -S'Cx,)  -  -  -(a  +  d-Z-.O) 

Sirre 


-  ncx^)  - 

:'(2a) 

-  J<D  -Ji)^  +  (x  -  b')  =71. 

*  (T-41-n) 

f(x) 

1 

(I-41-h) 

(“(x))" 

fCx^) 

•  ((Xj)  » 

1 

(I-41-c) 

''(Xj) 

=  (a 

+  1 
ca 

£211 

a-41-e) 

T)ie  first  order  solution  to  the  integral  is  given  by  (1-37)*  with  the  aid  of 
^T-40)  nnd  (I-41-a,  «•••  e)  as  follows 


^^n^flrsc 

order 


+  jkk„)C2T.) 
e _ 


(W‘3) 


Tc  obtain  tho  socon-i  order  terr-  (i.e.,  thnt  involvin';  f'(::)),  we  calculate 
t'le  I'ollowir.-;  quantities 


r'x) 


-n 


(I-43-a) 


(T_i3_bi 


; 
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V"(x) 

.  _  l!:p.A3ln  X  ^  3^^^^  ^ 

“1 

(I-43-c) 

f(Xj) 

■  -’'•''2^  """mT 
/U 

(1-4  3-d) 

<!>"(»l) 

■ 

(T-43-e) 

4"'(Xj) 

-  <h"'(x2)  -  0 

(T-45-f) 

Substitutins  (1-40)  and  (I-43-a,  ....  f)  into  (1-37)'  and  adding  thn  first 
order  solution  (1-4?)  the  Integral  to  second  order  is  given  approxinately  by 

(J  \  ^  ^  *  - 

'‘^n'^first  and  “  n  (a  +  jl;1>  )gi 
second  ^ 

order  terms 


l(a  +  i  >v 


1 

a  +  j>c'!i^„)V2J 


(T-/'4) 


'ore  pertiel  integration  steps  would  leed  to  ter-ns  with  increasingly 
higher  or:>r  variation  with  ^ ®nd  the  ratio  (^).  The  latter  rntio  approaches 
unity  and  }\, approaches  as  distance  fron  the  cable  decrenses.  lerce 

close  to  the  cable,  we  would  obtain  a  series  in  increasing  powers  of  the 
reciprocal  height  of  the  observation  point. 

eJ'”' 

'  e  note  that  the  factor  ^  is  present  in  this  solution,  this  Is  iincr- 
tant  in  supporting  the  argument  that  the  fields  near  the  cable  are  very  close 
to  those  that  would  he  obtained  fror  only  those  slots  nearest  to  ttc  ohservof'on 
point. 


II-l 


A')oendix  IT 

Petailg  on  Approaches  '.’sin".  Plnne-'  ave  Snectrni  Rooresentation  of  Fields 


In  C‘>e  orininal  nroiect,  cho  planc-uave  spectral  representation  of  each 
fiel'l  comonent  (See  lef.  1,  Sections  3  and  4  and  Appendices  II  and  III) 
or  a  valid  nporoxir.ation  to  it,  was  needed  for  two  purposes: 

(a)  To  determine  t''e  ground  reflected  fields,  based  on  the  assumption 
Qt  a  plane  wnve  incident  on  the  ground  surface 

(b)  in  order  to  validly  use  Che  Parber  scattering  progran  ("ef.  1, 
Section  71  which  assumes  plane  waves  incident  on  the  scatterer. 

The  theory  behind  the  transformation  between  position  space  (x,  y,  z 

or  (5,  d,  z  coordinates)  and  "plane-wave  spectrun"  space  (the  space  of  the  _ 

■ave  vector  '2’  ''hose  r,  y  and  z  coordinates  are  where  ■  ±  - 

is  given  in  'eferirce  1.  Appendi::  IF.  If  ^(x,  y,  z)  is  a  vector  function  of 
rectangular  coordinates  (x,  •/,  z),  then  from  F.q.  (II-3)  in  Appendix  II  of 
Sof.  I,  with  some  notalional  changes,  we  have  the  inverse  Tourler  transform 
of  ^(x,  y,  z)  in  terms  of  the  plane-wave  spectrum  for  upward  and  downward 

^  A  A 

nropagnting  wave?  an('  ''  (^)  respectively,  where  =  3„  j  cs 

follows: 


_  „2 
”  '  V 


i(x,  y.  z)  .  ••  •  1^  ;  Vja,) 


j'ih  Iz/V>  1 


,1  " 


1. 


Ve  wrill  nov  conbine  (TI-7-a,  b)  with  b)  or  ^2.10-a,  b)  to  fom 

the  integrands  of  the  expressions  <II-2)  where  the  vectors  V  and  V  in  (TI-2) 
are  the  electric  field  vector  ^  due  to  the  slot  at  ■  o  +  '3'  and  its  plane 
wave  spectru-i  ^  respectively. 


ft  -u*  r  J- 


2n  "  ®h  ^  ~ 


-(«  +  +  9’),  ^ 


f.  .,\  ■  f  |3  I  sin  •>'  ] 
Uj!'^  )  '  z'  cn 


+  - ^  '(JkP  cos  0  -  j!'jl  cos  0')'n  I  +  jhb’ ,  sin  *'] 

fJkR^)- 

4  .  f(jhp  'cos  -  j't^.cos  >i’)(-  1,1^1  ±  sin  o'  : 

— T  cos  -  j;'4(lcos*n*)f±  3 

4  - ^ — r  iCjbP  cos  •■)  -  jhAcos  i')(+  3  j'-.i'))|  (IT-l-r.) 

,  i  I"  f2ii  jU['’  -  3,,  o  cos  (t  -  c-1] 

(c?,(^)].Jo.  -  fip,  ^  ® 

'  ^  -  4n''  O  V, 


^ j  [(J'<0  sin  9  _  ji^.sin  i  JVi*  !  cos  oM 


(J':Ro)- 


fijho  si'n  9  -  sin  r’)'-  +  j':3' 


■  ’■■■<  •  'i 


1 


sin  n  -  j'Asin  ‘.'M-  b 
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=q 

47-2  JQ 


2t  jl-.;i’.  -  'j  cos  (;■;  -  *;)’ 


1  + 


-(.'.+  j!-'i\  )i(*  +  e-)  ,  ,  .-1 

e  -iV  f"i  ^7:::7-2  ii^.i  b' 


f-  i^J  j!^  +  - L_^  t;  i  (j;.)'  (b')2l 


^  [±  3  (jk)-  (h')2l| 


7jt-P  )-’ 


(!’■-?) 


The  quantities  fT^(3;^)]^,  whether  calculate:?  fro"  (TI-'— ci  or  (11-9) 
rust  be  exactly  equivalent.  Of  course,  fro"  a  practical  vie'.noirt,  the  nost 

^  A/ 

efficient  '/av  to  obtain  T  Is  to  cc-?jte  P  and  .  tron  '’’I-  -a,  b)  and 

2  *  V  y 

V 

t  i^n  confute  fron  (TI-j),  The  only  re'\son  one  nigiit  cornnte  ^ron 
is  for  checl;ing  purposes. 

Carryir."  out  the  integrations  around  the  cnbl-  slots  in  (:I-8-c,  b,  c) 
and  invoh.inn  (3.17-a,  b,  c)  and  the  enprohin-tiond’.l  t:e  ohtair 


■Y>>'= 


rafi  f>  -Jb 


n  cos  (o  -  hi 


1  *  13  ! 


oin 


(A  A  A  A  A  ^ 

^  h  ^  ‘’x  ^2  *  ^  ^3  •=•  ■  q  ^4  ■  ^  ^ 

V  V  y  V  v  ' 

z  1  i  ♦  * 


.'bore 


'3  I 

— jT-k'  sin 
ca 


f  T  T  _  7  O. ) 


^x  ■  TTO  '-z'  ‘  'ca 

c  «  — r  '(0  -  JO  cos  0  (-  i'  I  +  ji  b'l  t  b'  !:*  sin 
X  z  cr 


±  - 5“  (3b  (p  •  ^  cos 


+  - T  r3b*  (p  -^)  cos  P] 


'  7.  A, 

/  V  cos  •,> 
U*<)  ca 


=  TjTry  [(p  -A)  l'"^!  si*'  ^  ±  b*  cos  r^) 

=  — ^■  r(c  -  a)  sir  A  |2  I  ?  <Vd’)  ?  ■•'  cos 

'■  “  "- 

*  ^  ""  „  r ?h  *  sis 


»  -  r3h'  (P  -  ^)  sin  o} 

Mi-)-' 


(j'O  '"h 


sin  (  i  -  *„) 


'  ( j':)  I'rj  j  '’rJ  con  f  T  -  '*.)  +  '>’  -,1^  sin  (/>  -  (lip)] 

“  [(0  -  i!;i  cos  (■>  -  •.,)(*  +  jUV)  .  b'  sin  (i  -  <''^1 

=  - r -  '3b'  (r.  -  a)  cor,  -  *,>! 

X  \ 

=  - -  l’3b'  (o  -  cos  (-'•  -  i,)' 

'j'-)  |3  !  ■ 


2n  ,  -  f  j‘.J:  C-HiO’ 

JO  cc 


T,  =  I  iC-i-c 

r  J  ^ 

o 


-  /^"  +  iT  -  2  5^  COS.3' 


■\  serious  ’ttmn'  .os  iinie  tr  enrry  out  the  (’ouhir  intioyrntions  innicater 
n-'-n,  c),  ioc^-  ■':■  3pnror:*-'nte  ans’ytiCcsl  r.otliorts  nr.d  mnericnl  netbocis 

■'  rn  t’’.-"  -rtiio'  urs-i  t'>r  lnt''erotior  nro'.nr  the  r..n')lo,  ''Isc'.'ssad  in  ippcr.'H-; 
ec  1 1  on  1 -r ) . 


'  •**  i* fc'  •  ***>'*..r^*  *■**  *•**  J,*"* 


■-■.v 
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5^ince  t'le  are  independent  of  </,  and  tlie  remainder  o'  t'le  mnlitri'e 

of  t*ie  integrand  varies  slowly  «lth  the  intejratior  on  'i  appears  to  be 
very  easy.  Stationary  phase  nethods  were  considered,  bat  it  was  found  that 


the  stationary  phase  points,  piven  by 


sin  -  *,) 


(II-ll) 


do  not  exist  because  the  -"S  of  (II-ll)  exceeds  unity  over  n  '■.i5nlficant 
portion  of  the  intejration  renion.*  T!ie  other  approxination  rethod  discussed 
in  ‘ppentix  1,  Section  C,  applicable  in  cases  idiere  stationjry  niiese  points 
don't  exist,  was  also  considered.  Plnally,  the  f)  integrations  '..'ere  actually 
carried  out  using  the  very  nccurnte  nurerical  nethod  discussed  in  '.npendix 
I,  Section  ■?.,  :  os-3ver,  when  the  integration  on  9  wms  attcnptef  by  nur-erical 

nethods,  it  ••as  foend  that  coverage  of  the  reeion  near  the  cnhle  (the  eost 
irnortnnt  contribution  to  the  integral)  recuired  a  prohibitive  mount  of 
eonpiiter  tine  and  that  apnronch  was  ahandoned."  At  that  toint,  ?n  anprexinote 
approach  i-as  considered  wherein  the  field  conponents  are  given  hy  f’.'-a, 
b,  c)  only  within  a  region  that  contains  the  scatlorer  volune  and  vanish 
outside  that  region.  Thin  night  be  e  satisfactory  approach  becaiso  the  scetterer 
in  free  space  only  responds  to  the  fields  irpinging  tipon  it  and  the  scatterin'* 
is  independent  o'  the  liel'is  in  the  surrounding  raglons.  Of  course,  the  plana 
wnve  sprectru'i  obtained  hy  tais  techninTic  \oul  '  he  <li  '  'ernnt  fron  f’at  ohta'ncd 
by  integratir.f’  (I''-!'))  over  oil  space  (e.g.,  ther';  unulc  be  si--.ificant 

*  11. b 

Ta  order  for  r.  stationary  nhase  aoiat  to  o?'ist,  o  would  hnvo  to  excee''  -rr—  ■ 

•h 

i^ince  me  lar'^ont  aosslhle  value  of  d,  i;;  I,  that  would  on’-  aanlv  "or  nointr 

h  ’  .... 

frr  outside  the  cable. 
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contribution-;  fro--,  a  larf.er  portion  of  ^-space).  However,  when  the  inverse 
Fourier  transfonntion  is  nerforoed,  the  recovered  field  conponents  should 
be  approximately  the  same  as  those  obtained  throii<'.h  integration  of  fIT-10) 
over  all  position  space. 

The  numerical  inplementation  of  this  nethod  is  very  sinole.  It  anounts 


to  replacin3  the  double  Intenration 


f”  f2it 

dPP  •' 
•'O  ^G 


dA  in  (II-IO)  by 


O  i  ^ 

fo  +  2 

^^o-f 


To** 

j  ^  ''q)  is  the  center  of  the  inteoration  volume,  ip 

0  2 

t.ne  width  of  that  Integration  voU'ne  in  p-space  and  Ao  f-.e  wlP.th  in  A-snnee, 
'Tie  above  apnroech  was  not  actually  Innlemcnted,  Although  it  would 
have  been  ideal  to  attain  the  plane-wave  spectral  field  components  by  actunll; 
c-rryin3  o-jt  tin  intesrations  in  (Il-^-a,  b,  c)  there  was  Insufficient  tine 
ani  rosourcas  to  do  Justice  to  that  aporoach.  It  was  then  decide-'  to  adopt 
an  mnronch  that  was  also  used  in  the  original  project  and  -diich  circimvonts 
the  need  for  the  double  intesratioa  in  (IT-IO).  That  approach  is  based  on 
the  relationship  (III-l)  in  leference  1  derived  in  Appendix  Tllof  'lef.  1 
nn-l  r“ponted  below  with  snail  n-itational  modifications. 

(1.-  i'.p  f2if  fS  I'  ,  -I  jlt-T  cos 


■  tiore  thc’int-^r^cion  over  the  -jo'  ericnl  rinr^ies  (0^,  covers  nil  pos.'^iMe 
r^irectiom  the  ^-vector, 

i'i  idcnte*!  to  tho  In  this  sturt;*  hy  firPt  notin'* 


■«  III  Ilf  M.  l» 


; 


11-10 


f  dO,  sin  €,  .  t 

Jq  ■'  Jo  ^  ^  Jo  “  ^ 


(IT-13-b) 


Usln:;  (II-13f-a,  b)  in  fll-12>,  i.e  obtain  the  latter  in  the  rorst  • 

1!^.; 


i’^r  1  f“^ 


ni-u) 


or  equivalently,  In  order  to  eeparste  out  "upward"  and  "downward"  oropa-atlnj 
waves 


-:i:2 

(?“ 


,  [In  /I  Ih3.  (p  -  (l)  cos  (*<  -  h.  ) 

^lo 

,  jk|3  |(z  -  b') 

[e  *  ■  (1  +  ji/n  +  i'.;)) 

-Jk|3j(a  -  •)•)  _  , 

^  ^  -  JW  J-:(n  «  ^  .;}) 


'.’here 


f  I'^J  »  \  <i  cos  -h  I  sin  -I,)  -  $  ]\\ 

r  ■  i  fix  '•  i  ■  i  I 'si  “  ‘’h  ‘='’®  ''3  +  X  ■,)  -  i 


■’,  ' 
's' 


The  orocedure  based  on  the  use  of  (11-16)  hefins  by  nsnunln",  a  "locnl 
plana  wave  spectru-i."  "y  "local,"  '<e  necn  that  the  planc-'.'avc  spectral  field 
co-iponents,  which  are  indeperdonr  of  position  coordinates,  are  appraxl'iated 
at  a  point  in  position  space  hy  the  actual  field  coiponents  (2.1'-n,  a,  c) 
cultlplied  by  a  particular  factor.  That  factor  is 

£-(^)  -  ^  ^  +  Jh(R(C)  +  |-.J(0))j  c-J  c  *  (TI-15) 

where  2(0)  •  distance  between  observation  point  and  the  nearest  cable  slot  ■ 
/(p  -  e  (z  -  b')^  and  j'(ii)  ■  position  of  nearest  calle  slot  «  cos  d;) 
+  sin  ■b)  +  jh' . 


■;-y 


I  . 


c 


i_ 


’V-'\ 


.^A.TjvWLv\i 


II-ll 


T!i“  lo^if  hehlntf  thin  -'.pproach  1*  axplalned  below.  The  field  cortponents 
fron  an  individual  slot,  as  3iven  by  (2.3-a,  b,  c),  all  contain  the  factor 

idiere  ?  Is  the  distance  between  that  slot  and  the  observation  point. 
i':R 

■nie  quantity  o-' ' "  is  then  expressed  as  a  superposition  of  plano-'.'ave  spectra 

throinh  'q.  (71-14).  In  the  oriqinal  project,  this  was  done  for  each  slot 

and  then  the  Integration  over  the  slots  was  carried  out  for  the  "local"  spectral 

ivq 

field,  w!  ere  the  factor  e"'  "  v/as  replaced  by 

(  i  j’***  - 

^  [1  +  J1--  (”  h 

The  integration  was  then  accomplished  hy  approxinate  methods  (Ref.  1,  Sections 
3  and  4  in  Appendix  III). 

In  the  present  project, the  intecration  over  the  cable  was  done  for  the 

actual  fields  as  functions  of  position  coordinates,  not  for  their  plane-wave 

spectra.  Therefore,  havin-;  acconpltshed  that  integration  and  obtaine<)  thn 

field  components  fron  the  entire  cable  at  a  point  in  space,  we  would  li!:e 

to  irolenent  the  "local  plane  wave  r.pectruis"  model  based  on  these  results 

(i.e.,  based  on  Tes.  (2.’9-n,  b,  c)  rather  then  (2.C-a,  b,  c)).  Rowever, 

the  deterrent  to  that  is  the  fact  that  the  factors  e-  '  have  disaoneered 

pji-'  ’ 

in  the  process  of  intepratinq  around  the  cable,  'the  factors  -= -  in  f2.j- 

■  •  -  .,n  ■ 

a,  b,  c)  or  (2.10-a,  h,  c)  have  been  replaced  by  in  (2.1?-a,  I),  c>). 

If  is  evaluated  by  the  partial  integration  method  indicated  in  Section 

d-f'  of  Appendix  I,  the  factor  appears  in  T^.  lowevcr,  since  the 

more  accurate  numerical  intep,rction  technique  discussed  in  Section  I-C  of 

Aopendi:;  I  "as  nctiiallv  used  tc  obtain  I  ,  the  factor  e-*  '  '  '  is  not  explicitly 

n 

indicated  but  wa  know*  that  I  could  be  written  In  the  fom 

n 


■"  ”bo  results  using  the  partial  '.ntcgratioa  method  were  compared  with  tboaa 
from  the  numericnl  technique  and  the  apreemont  was  nearly  perfect. 


.v.v.-* , 


•  •>  ■_» 


n-12 


-  1 

n 


(II-IS) 


A  field  component  as  jivcn  by  (2.19-a),  (2.19-b)  or  (2.19-c>  could  be 
written  (throush  (11-16)  and  (TI-14)  in  the  forn 
3  /  3 


5(r)  -  I*:  I  -  fjElj  e-*'-'^ 

^  n-1  "  "  ^n-l  " 


-R(0) 


,2ir  ,1 

h  'G 


-h « 


%.  1) 


(IT-17) 


where 


^(a.  i) 


1  I  V  -n-t'COr 
^  I  r  e  !l  ^  iurm 

n«l  *- 


•■  a.;(o»] 


[•i  r/r 
l4it  '* 


) 


Tbe  field  component  -(J,  g),  a  f'mction  of  both  ^  and  th?  observation 
point  vector  g  ■  (p,  P,  z),  is  interpreted  as  the  "locel  plene  wave  spectrin" 
of  ’(r)  and  is  obteined  by  rultipTylnp  the  nctuel  r/f)  by  tbe  factor 
^  (I  +  jV.C'TOI  +  i.J(0))]e  e  as  Indicated  in  (TT-15). 

This  method  was  used  and  a  niinber  of  corpatrtiona  ■•ere  perfomed  with 
it  including  prod'iccion  runs  for  scattered  fields.  Tae  results  were  Jud:;cc 
to  be  unsatisfr.ctory  and  eventually  i-t  ••as  abandoned  nn-l  replaced  byt'ie  -ethod 
described  in  Section  4  in  the  msin  body  of  this  report. 

One  of  the  major  deficiencies  of  the  nethod  is  as 'foUo'JE.  "hile  it 


produces  a  spectrum  of  "plane  l•aves,"  i.e.,  waves  'rttb  the  plane  wave  propepo- 
tion  factor  e  ,  there  wave?  are  not  necessarily  trrnsverve.  The  oethod 


described  in  "ectioa  4,  or  whic'-  our  n'jrerical  r^nsvlts  arc  based,  do"s  produce 
an  appro::lnntion  to  the  true  ^leld  ct  tbe  scatterer  thet  ren^  '  s  TT! 

•vave. 


.’’Vv V,  .rVL-/>.T.v 


fn  the  table  below,  the  nAr^^neters  inputted  into  the  conputer  nro^ran 
are  list?'*,  v/ith  the  algebraic  synbol  uned  ''or  each  vorirble  li  analysis, 
the  Fortran  nme  of  the  variable,  its  dc^’inition  and  the  values  assigned 
ro  it  in.  the  pro^ran.  This  Ttaterial  ir?  ti!:cn  fron  ^lef.  1,  Section  11,  Papes 
1)-*?  thrnurh  11-6.  The  oripinnl  raterinl  lias  been  changed  to  reflect  the 
if  ice  r.  ions  nnde  in  the  rev'  proyren. 


Paranetera  associated  with  the  cable 

Fort  nn  My.ehrnic 

lane  Fvnhol  DcTinition  ^>lue 

''adlus  of  cable  24 

configuration  in  ^nonroxination 
r’eters  based  on  circiin- 

forenco  of  151 
nsters) 


\  a  inner  radius  of  .00470 

coaxial  caole 
r.eters 


b  Outer  radius  nf  .0127 

coaxial  cablo- 
notors 


'.'■n-'t Slot  length  (alon’  .003 

cable);  s.t?c  for  all 

slots-meters 


Slot  width  (around 
cable);  sane  'or  all 
slots-faeters 


0155 


cable  naterlr.l 
n*»os/ii«ter 


Azinuthal  angle 
along  the  cable 
of  Slot  ’1.  (The 
slot  nearest  the 
pou-er  source)  in 
eround  frane 


An-'le  of  slot 
center  around 
periphery  o^ 


/.ttenuation  or 
Ti:  node-nepers/ 
netcr 


’adlo  frequency 
in  ’iertz 


57(10"} 
57  :!!:* 


Voltage  bet veer 
Inner  and  outer 
conductors  of 
coaxial  cable- 
iisec  in  arplitudn 
of  TC  -loce-volts. 
?arar>ctcrs  associated 
with  the  ".round 


Per-.ittiv’ ty  of 
ground-*^"ra''s/ 
'■etcr 


iroun"  con'  i’c- 
Li  Vi  ty-'hf'os/'aeter 
^nrrrreters  assoc¬ 
iated  with  the 
reorpt  rv 


a-coordinate  o’ 
•antenna--ef'r''. 


'  "r.’."nnn  at 
center  of  erreeJ-r 
conti-iThtior } 


y-coordinato  of 
antenna-fneters 


(same  renar’cs 
as  above) 


T!!B 

Spherical  polar 

angle  of  long 

axis  of  prolate 

spheroidal 

scatterer 

(in  projrd-fr.ine) 

Varied: 

0  i'or  "upright 
nan"  cases; 

y  for  "crav;linp 
r.an"  cases 

pTjr. 

'b 

Azimithal  angle 
of  long  axis  of 
scatterer  (in 
groutid-frcTc) 

Var  iod 

C 

« 

Pernittivity  of 
scatterer  rater ial- 
farads/netcr 

A'l.TCPClO"'-) 

4'i.2E^ 

0 

hosed  on 
'eferenco  h.ic 

SI'T'AS 

Conductivif;  of 
scatterer  rateritl 
"ihos/neter 

O.i''-’ 

ho.sed  on 
■'oCerence  47 

)C 

■•s 

1  2  ^ 

W  +  y*;  «  rn^’lel  Vari5(. 

(ccr'.wte<' 

v-.rls)!? 

coordinate  of 
scatterer  center 
in  ground  frane 

P’IS 

(co^pute<^ 

vnrinbls) 

'*5 

S 

-1  f’'s1 

tan  — '  •  azlButh  Vnrieri 

angle  of  scatterer 
center  in  ground  fra-ie 

r.s 

o 

■■5 

Radius  of  nrolatc 
spheroidal 
scattercr-net  er s 

:.r 


i  of  proljt? 
>;pi'9roi''ol 
sect  ter or-r  (•tern 
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